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ABSTRACT 


The  present  investigation  is  a  study  of  the  modeling  of  radial 
flow  gas  turbine  fluid  dynamics  and  heat  transfer  behavior.  A  test 
facility  has  been  constructed  for  the  modeling  of  turbine  flows. 
Several  unconventional  experimental  techniques  have  been  developed 
for  use  in  the  facility. 

A  simple  model  of  a  radial  flow  gas  turbine  has  been  analyzed 
and  tested.  The  fluid  dynamics  and  heat  transfer  behavior  of  a  single 
radially  inward  flow  of  air  over  a  shrouded  rotating  disk  has  been 
experimentally  determined.  The  fluid  dynamics  behavior  is  presented 
as  a  series  of  graphs  detailing  the  velocity  profiles.  The  heat  trans¬ 
fer  is  correlated  over  the  operating  range  of  radial  flow  gas  turbines 
by  a  simple  algebraic  relation  between  the  governing  nondimensional 
parameters  .  An  analysis  is  presented  that  does  not  adequately  pre¬ 
dict  the  heat  transfer.  The  reasons  for  the  discrepancy  between 
analysis  and  experiment  are  discussed. 
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NOMENCLATURE 


English  Letter  Symbols: 

Cp  fluid  specific  heat ,  BTU/lb  F 

2 

g  gravitational  acceleration,  ft/sec 

2 

h  average  heat  transfer  coefficient),  BTU/hr  ft  F 

J  (x)  nth  order  Bessel  function  of  first)  kind  with  argument  x 

K  correlation  parameter  defined  by  Eq.  [17] 

heat  meter  calibration  factor,  (feTIJ/hr 

mv 

k  thermal  conductivity,  BTU/(hr  ft^  F/ft) 

9 

p  fluid  pressure,  #/in 

2 

q"  average  heat  flux,  BTU/hr  ft 

r  radial  coordinate,  ft 

rQ  disk  radius,  ft 

T  temperature ,  R 

t  temperature ,  F 

u  fluid  radial  velocity  component,  ft/sec 

V  fluid  resultant  velocity,  ft/sec 

v  fluid  tangential  velocity  component,  ft/sec 

W  mass  flow  rate,  lb/hr 

w  fluid  axial  velocity  component,  ft/sec 

z  axial  coordinate,  ft 

zQ  disk-to-shroud  spacing,  ft 


xiii 


Greek  Letter  Symbols: 


2 

a  thermal  diffusivity,  ft  /sec 

a  m  eigenvalues,  zeros  of J^(x) 

g  Reynolds  number  ratio; 

air  coefficient  of  thermal  expansion,  1/F 

6  hydrodynamic  boundary  layer  thickness, 

A  thermal  boundary  layer  thickness,  ft 

O 

eddy  diffusivity  for  heat,  ft  /sec 
H 

o 

eM  eddy  diffusivity  for  momentum,  ft  /sec 

n  dummy  variable  of  integration 

6  tangential  coordinate,  radians 

y  fluid  dynamic  viscosity,  lb/hr  ft 

o 

v  fluid  kinematic  viscosity,  ftv  sec 

p  fluid  density,  lb/ft^' 

2 

t  shear  stress,  #/ft 

o  disk  angular  velocity,  rad/sec 

Nondimensional  Parameters: 


NGr 

9  3  ro 

p  Qrashof  number 

Z 

u 

NNu 

h  ro 

average  Nusselt  number 

k 

nr,f 

W 

ZQ  y 

mass  flow  Reynolds  number 

nr,m 

2 

r 

o  u  p 

,  machine  Reynolds  number 

U 

xiv 


N. 


St  /  W  ,  ^ 

2  tt  r  z  °P 
o  o 


,  average  Stanton  number 


N 


St,  1 


W 


, - -)  cD 

2 *  ro  zo 


,  local  Stanton  number 


Subscripts: 


o 

z 


nondimensional  radial  coordinate 


nondimensional  axial  coordinate 


]o_  ,  nondimensional  disk-to- shroud  spacing 


Lo 
2  it  N 


N 


,  ratio  of  machine  Reynolds  number  to 


R,F  mass  flow  Reynolds  number 


,  nondimensional  hydrodynamic  boundary 
layer  thickness 


— -  ,  nondimensional  thermal  boundary  layer 


thickness 


_A _  ,  boundary  layer  thickness  ratio 


r  denotes  value  in  the  radial  direction 

R  denotes  relative,  evaluated  with  respect  to  the 

rotating  disk 

6  denotes  value  at  the  edge  of  the  hydrodynamic 

boundary  layer 


xv 


m 


denotes  fluid  inlet  value 


1  denotes  local  value  as  opposed  to  the  average  value 

Superscripts: 

*  denotes  a  nondimensional  quantity 
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#  denotes  pounds  force 

lb  denotes  pounds  mass 
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I  INTRODUCTION 


In  recent  years  there  have  been  increasing  applications  of  radial 
flow  gas  turbines  as  prime  mover  and  auxiliary  power  sources  in  the 
aircraft  and  automotive  industries.  High  performance  levels  in  these 
units  are  attained  with  high  turbine  inlet  temperatures.  These  high 
temperatures  necessitate  the  use  of  high  temperature  alloys  in  the 
turbine  wheels.  A  proposed  method  of  reducing  the  metallurgical 
requirements  is  to  insulate  the  turbine  wheel  from  the  hot  gas  stream 
by  a  veil-flow  blanket  of  cold  air,  as  shown  in  the  following  figure: 


A  heat  transfer  investigation  of  the  veil-flow  method  has  been 
undertaken  jointly  by  the  AiResearch  Manufacturing  Division  of  the 
Garrett  Corporation,  Los  Angeles,  California,  and  by  Stanford  Univer¬ 
sity.  AiResearch  has  conducted  a  prototype  study  of  a  veil-cooled 
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radial  flow  turbine  (16)  .  The  present  investigation  is  concerned 
with  model  studies  of  radial  flow  turbines. 

The  present  model  consists  of  a  shrouded  rotating  disk  to  which 
vanes  may  be  added.  The  main  air  stream  flows  radially  inward  from 
the  disk  periphery.  The  veil-flow  is  introduced  next  to  the  disk  at  the 
periphery,  and  also  flows  radially  inward.  The  disk  and  veil-flow  are 
heated  relative  to  the  colder  main  stream.  The  heat  transfer  results 
obtained  from  the  model  tests  will  be  extended  to  radial  flow  turbines 
where  the  main  flow  is  hot  relative  to  the  colder  disk  and  veil-flow. 

Because  of  the  complex  flow  and  heat  transfer  behavior  expected 
in  the  veil-flow  studies,  the  initial  investigation  is  concerned  with  a 
single  air  stream  flowing  radially  inward  over  a  shrouded  unvaned 
rotating  disk.  The  fluid  dynamics  and  heat  transfer  behavior  for  this 
model  will  serve  as  a  datum  for  the  later  investigations. 

The  study  of  a  shrouded  rotating  disk  is  also  of  general  interest. 
Many  heat  transfer  studies  have  been  conducted  on  unshrouded  disks, 
and  several  on  shrouded  disks  with  radial  outflow.  The  present 
investigation  is  another  phase  of  the  rotating  disk  problem. 

The  present  investigation  has  as  its  specific  objectives  to: 

1.  Design,  construct,  instrument,  and  evaluate  a  test  facility 
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Numbers  in  parentheses  refer  to  references. 
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for  determining  the  fluid  dynamics  and  heat  transfer  behavior 
of  a  rotating  disk  for  various  flow  configurations. 

2 .  Analytically  and  experimentally  obtain  the  fluid  dynamics 
and  heat  transfer  behavior  for  a  single  radial  inflow  of  air 
over  a  shrouded  unvaned  rotating  disk. 

The  present  investigation  has  satisfied  the  first  objective  and, 
to  some  extent,  the  second.  A  test  facility  has  been  constructed 
that  allows  testing  of  various  rotating  disk  flow  configurations. 

Several  unconventional  experimental  techniques  have  been  developed. 
The  disk  is  heated  by  a  controlled  infrared  radiation  source.  Heat 
meters  are  used  to  measure  the  heat  flow  through  the  disk.  A  novel 
method  of  transferring  the  thermocouple  and  heat  meter  electrical  out¬ 
puts  from  the  rotating  disk  to  stationary  recording  instruments  has 
been  developed. 

The  rotating  disk  fluid  dynamics  and  heat  transfer  behavior  has 
been  experimentally  determined  over  the  operating  range  of  radial  flow 
gas  turbines.  The  flow  pattern  is  unexpectedly  complex,  showing  a 
strong  influence  of  the  presence  of  the  disk  and  the  shroud.  The  heat 
transfer  results  are  correlated  by  a  simple  algebraic  relation  between 
the  governing  nondimensional  parameters.  An  analytical  attempt  to 
predict  the  heat  transfer  is  also  presented.  The  discrepancies  between 
analysis  and  experiment  are  discussed. 
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II  SUMMARY  OF  WORK  OF  PREVIOUS  INVESTIGATORS 


A  survey  of  the  available  literature  pertinent  to  the  model  study 
was  undertaken  to  help  satisfy  the  project  objectives.  Considerable 
information  exists  for  heat  transfer  from  unshrouded  disks,  but  very 
little  for  shrouded  disks.  There  is  also  some  information  on  heat 
transfer  in  veil-flows  over  flat  plates  which  may  possibly  be  extended 
to  veil-flows  over  rotating  disks.  A  comparison  of  the  analytical  heat 
transfer  results  in  the  literature  has  proved  helpful  in  evaluating  the 
various  analytical  methods. 

A.  Heat  Transfer  from  Unshrouded  Rotating  Disks 

Numerous  studies  of  the  heat  transfer  from  an  isothermal 
unshrouded  rotating  disk  to  an  infinite  medium  have  been  conducted. 
The  unshrouded  disk  geometry  and  flow  pattern  is  shown  in  Fig.  1, 
and  the  heat  transfer  results  are  summarized  in  Fig.  2.  The  experi¬ 
mental  laminar  heat  transfer  results  of  Cobb  and  Saunders  (2)  and 
Kreith,  Taylor,  and  Chong  (12)  are  in  good  agreement  with  the  ana¬ 
lytical  results  of  Dorfman  (4),  Millsaps  and  Pohlhausen  (13),  and 
Wagner  (22).  The  experimental  results  of  (2)  are  obtained  directly 
from  a  heated  disk,  while  those  of  (12)  are  established  by  applying 
the  heat  and  mass  transfer  analogy  to  the  mass  transfer  from  a 
napthalene  disk.  The  agreement  between  analysis  and  experiment 
in  the  laminar  regime  establishes  confidence  in  the  experimental 
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results  in  the  turbulent  regime. 


As  shown  in  Fig.  2,  there  is  considerable  disagreement  in  the 

turbulent  regime  between  the  analytical  attempts  of  Davies  (3), 

Dorfman  (4),  and  Kreith  and  Taylor  (11),  and  the  experimental  results 

of  Cobb  and  Saunders  (2),  Kreith,  Taylor,  and  Chong  (12),  and 

Zaloudek  (23)  .  It  should  be  noted  that  all  of  the  experimental 

results  contain  the  effect  of  a  laminar  core  at  the  disk  center;  the 

flow  does  not  become  turbulent  until  a  local  Reynolds  number, 

2  5 

(r  w  p/p  ),  of  2  x  10  is  reached.  At  a  machine  Reynolds  number  of 

r 

10  ,  the  flow  is  turbulent  over  80%  of  the  disk  surface.  As  the  ana¬ 
lytical  solutions  are  for  turbulent  flow  over  the  entire  disk,  the 
experimental  results  should  approach  the  analyses  as  ^  increases. 

The  result  of  Davies  (3)  was  obtained  from  an  integral  analysis 
using  a  very  approximate  evaluation  of  the  heat  and  momentum  dif- 
fusivities.  The  analysis  of  Dorfman  (4)  is  based  on  the  von  Karman 
analogy,  where  the  universal  velocity  profile  for  turbulent  boundary 
layer  flows  is  assumed  and  the  turbulent  diffusivities  for  heat  and 
momentum  are  equated.  Kreith  and  Taylor  (11)  present  two  analyses. 
Their  integral  analysis  result  is  obtained  from  the  integrated  energy 
equation.  The  hydrodynamic  and  thermal  boundary  layer  thicknesses 
are  assumed  equal.  The  analytical  hydrodynamic  solution  of  von 
Karman  (8)  is  employed.  In  their  Reynolds  analogy  result,  the  sum 
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of  the  laminar  and  turbulent  diffusivities  for  heat  and  momentum  is 

assumed  equal  over  the  entire  boundary  layer.  The  Stanton  number 

is  thus  equal  to  one-half  the  friction  factor.  All  of  the  analytical 

0.8 

solutions  are  of  the  same  form,  N^u  =  C  ^  ,  with  the  value  of 

C  dependent  on  the  analysis  employed. 

The  data  of  Cobb  and  Saunders  (2)  and  Kreith,  Taylor,  and 
Chong  (12)  are  in  good  agreement.  The  results  of  Zaloudek  (2  3)  are 
obtained  for  a  constant  heat  flux  surface  condition.  On  the  basis  of 
turbulent  flat  plate  heat  transfer  experience,  the  heat  transfer  coeffi¬ 
cients  for  constant  heat  flux  would  be  expected  to  be  about  5  per  cent 
higher  than  those  for  constant  surface  temperature,  instead  of  15  to 
20  per  cent  higher  as  reported. 

The  analysis  of  Dorfman  (4)  and  the  Reynolds  analogy  result  of 
Kreith  and  Taylor  (11)  compare  best  with  the  experimental  results. 

The  analysis  of  (4)  is  the  most  thorough  and  the  assumptions  made 
concerning  the  velocity  profiles  and  the  diffusivities  are  the  most 
plausible.  In  both  of  the  results  of  (11),  the  hydrodynamic  and  thermal 
boundary  layer  thicknesses  are  assumed  equal.  This  is  evidently  a 
poor  assumption  for  the  integral  method  approach,  and  casts  doubts 
on  the  accuracy  of  either  method  when  applied  to  the  shrouded  disk 
problem.  The  relation  of  Davies  (2)  is  in  disagreement  due  to  an 
inaccurate  representation  of  the  eddy  diffusivity. 
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The  best  analytical  approach  to  use  in  solving  shrouded  rotating 
disk  problems  appears  to  be  the  von  Karman  analogy  approach.  How¬ 
ever,  this  is  the  most  complicated  mathematically.  If  the  proper  rela¬ 
tion  between  the  hydrodynamic  and  thermal  boundary  layer  thicknesses 
could  be  found  to  make  the  results  of  the  two  methods  compatible,  both 
results  of  (11)  would  probably  lie  closer  to  the  data.  Unfortunately, 
this  approach  is  mathematically  impossible  for  the  l/7th  power  profiles 
chosen;  however,  it  is  possible  for  cases  involving  a  shrouded  disk. 
The  general  conclusions  from  this  comparison  are  that  if  reasonable 
velocity  and  temperatures  profiles  are  chosen,  all  methods  determine 
the  Reynolds  number  dependency  to  the  correct  power  and  evaluate  the 
constant  to  within  an  order  of  magnitude.  For  the  shrouded  disk  prob¬ 
lem,  a  combination  of  the  Reynolds  analogy  and  integral  method 
approaches  was  employed. 

B .  Heat  Transfer  from  Shrouded  Rotating  Disks; 

Considerably  less  information  is  available  for  shrouded  rotating 
disks.  Soo  (20)  has  determined  analytically  the  laminar  flow  pattern 
for  radial  in-  and  outflow,  and  Soo,  Besant,  and  Sarafa  (21)  have  ana¬ 
lyzed  the  corresponding  heat  transfer.  Kreith,  Doughman,  and 
Kozlowski  (10)  have  determined  the  shrouded  disk  radial  outflow  heat 
transfer  experimentally  using  the  analogy  between  heat  and  mass 
transfer.  Smith  (IS)  has  experimentally  evaluated  the  heat  transfer 
from  a  shrouded  disk  both  for  turbulent  radial  in-  and  outflow,  but 
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the  accuracy  of  the  results  is  questionable  due  to  the  experimental 
technique  employed.  A  comparison  of  the  present  radial  outflow 
results  with  those  of  (10)  and  (19)  is  presented  in  Fig.  34. 

C .  Heat  Transfer  in  Veil-flows 

The  AiResearch  investigation  of  a  veil-cooled  radial  flow 
turbine  is  summarized  in  (16).  Temperatures  on  the  turbine  wheel 
were  measured,  and  the  veil-cooling  flow  temperature  effectiveness 
was  established.  Sizable  turbine  wheel  temperature  reductions  were 
obtained  for  low  veil-cooling  flows  relative  to  the  main  flow.  It  is 
hoped  that  the  veil-flow  model  studies  will  provide  a  more  complete 
understanding  of  these  results. 

There  have  been  several  experimental  heat  transfer  investiga¬ 
tions  of  turbulent  veil-flows  over  flat  plates.  Seban  (18)  reports 
temperature  decays  for  a  heated  veil-flow  bounded  by  an  adiabatic 
wall  and  a  cold  free  stream .  Heat  transfer  coefficients  for  the  heat 
transfer  between  a  constant  heat  flux  wall  and  a  veil-flow  for  various 
veil-flow-to-main- stream  velocity  ratios  and  for  equal  initial  veil 
and  main  stream  temperatures  are  also  presented.  Hartnett,  Birkebak, 
and  Eckert  (7)  summarize  and  compare  the  temperature  decay  results 
of  several  investigators.  Heat  transfer  coefficients  for  the  heat 
transfer  between  a  constant  heat  flux  wall  and  a  veil-flow  with 
various  ratios  of  veil-to-main- stream  temperature  for  a  fixed  velocity 
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are  also  presented.  These  results  will  be  of  value  in  helping  inter 
pret  the  experimental  veil-flow  results  and  in  formulating  analyses 
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FIGURE  1 


ROTATING  DISK  COORDINATE  SYSTEM 
AND  GEOMETRIES 


(a)  ROTATING  DISK  COORDINATE  SYSTEM 

(b)  COORDINATE  SYSTEM  AND  FLOW  PATTERN 
FOR  AN  UNSHROUDED  ROTATING  DISK- 
WITH  RADIAL  OUTFLOW 

(c)  COORDINATE  SYSTEM  AND  FLOW  PATTERN 
FOR  A  SHROUDED  ROTATING  DISK 

WITH  RADIAL  INFLOW 
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FIGURE  2 


COMPARISON  OF  HEAT  TRANSFER  RESULTS 
FOR  UNSHROUDED  ROTATING  DISK  HEAT 
TRANSFER 

Average  Nusselt  number  NNu  plotted  as  a  function  of  machine 
Reynolds  number  N^  .  The  numbers  in  parentheses  corre¬ 
spond  to  the  references  from  which  the  heat  transfer  results 
are  taken . 

Laminar  Regime: 

(2)  Cobb  and  Saunders,  experimental 

(4)  Dorfman,  analytical 

(12)  Kreith,  Taylor,  and  Chong,  experimental 

(13)  Millsaps  and  Pohlhausen,  analytical 
(22)  Wagner,  analytical 

Turbulent  Regime: 

(2)  Cobb  and  Saunders,  experimental 

(3)  Davies,  analytical 

(4)  Dorfman,  analytical 

(11)  Kreith  and  Taylor,  analytical 

(12)  Kreith,  Taylor,  and  Chong,  experimental 
(2  3)  Zaloudek,  experimental 
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FIG. 2 


Ill  GOVERNING  NONDIMENSIONAL  PARAMETERS 

The  ultimate  aim  of  the  present  project  is  to  model  radial  flow 
gas  turbine  heat  transfer  behavior.  To  correlate  the  model  and  proto¬ 
type  behavior,  the  correspondence  between  the  governing  nondimen- 
sional  parameters  of  each  field  must  be  established.  The  gas  turbine 
parameters  are  established  in  the  literature.  The  rotating  disk  param¬ 
eters  must  be  determined  from  the  governing  differential  equations 
and  boundary  conditions. 

The  rotating  disk  coordinate  system  is  shown  in  Fig.  1.  The 
differential  equations  and  boundary  conditions  are  written  for  a 
stationary  control  volume  located  in  the  fluid.  The  disk-to- shroud 
spacing  is  assumed  small,  so  that  the  fluid  pressure  is  a  function  of 
radius  only.  The  heat  transfer  is  assumed  to  be  in  the  axial  direction 
only.  The  fluid  entering  at  the  disk  periphery  is  assumed  to  have  one¬ 
dimensional  velocity  and  temperature  profiles.  The  fluid  and  disk 
tangential  velocities  are  equal  at  the  disk  periphery. 

Continuity  equation: 

8  u  +  u  -i-  8  w  =o  [l] 

8  r  r  8  z 

Momentum  equation  in  the  radial  direction: 

2  2 

u  AGL  +  w  A_JG  -  _ u_  _  JJc_  dp 

8  r  8  z  r  y8z2  p  dr 
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Momentum  equation  in  the  tangential  direction: 
3  v 


3  r 


3  v  .  uv 
+  w  -  +  - 

3  z  r 


32  v. 


[3] 


3  z 


Energy  equation: 


3  t 
3r 


+  w 


3  t 
3  z 


32  t 


3  z 


[4] 


The  boundary  conditions  corresponding  to  bq>th  radial  flow  gas 
turbines  and  the  present  model  are: 


Atr=r^ 


u(rQ,  z) 


v(rQ,  z) 
w(rQ<  z) 
t(rn,  z) 


W 


2  tt  ^  zo  P 


in 


At  z  =  0  •' 


u(r ,  0) 
v(r,  0) 
w(r ,  0) 
t(r,  0) 


0 

r  w 
0 

4w 


The  variables  p,  u,  v,  w,  and  t  and  the  coordinates  r  and  z 
are  normalized  as  follows: 


u(rQ,  z) 


v(rQ,  z) 


W 


^  v  ro  zoP 


O  U) 
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u(r0,  z) 
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The  normalized  differential  equations  [  1  ]  through  [  4  ]  become: 


Continuity  equation: 

3  u*  +  u*  ^  /  rn  \  3  w*  =  0  [53 

3  r*  r*  \zoJ  3  z* 


Momentum  equation  in  the  radial  direction: 


3  u* 
3  r* 


+ 


ro 


w” 


3  U* 
3  z* 


W 


W 


2*  t0  z  p 


2  IT 


dp* 

dr* 


[6] 


Momentum  equation  in  the  tangential  direction: 


3  v* 
3  r* 


w* 


3  V* 
3  z* 


u*  V* 


2  TT 


[7] 
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Energy  equation: 


There  are  seen  to  be  four  groups  that  aovern  the  fluid  dynamics  and 


heat  transfer: 


nondimensional  geometry  for  the 
rotating  disk 


W 


z 


o 


u 


A 


nondimensional  flow  rate;  a  flow 
Reynolds  number 


[ _ W _ ) 

\2  "  ro  zoP  / 


=  6 


ratio  of  disk  peripheral  velocity 
to  the  entering  fluid  radial  velocity 


By  defining  the  conventional  disk, 
or  machine,  Reynolds  number  as: 


N 


R,  M 


2 

rO  <JL>  P 

u 


the  velocity  ratio  can  be  written  as: 


and 


N 


Pr 


nR  M  4 
NR,F 


ft 


the  fluid  Prandtl  number 


The  correspondence  between  the  gas  turbine  and  the  rotating 
disk  parameters  can  be  established.  Due  to  the  nature  of  the  model, 
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there  are  no  disk  parameters  relating  to  power.  This  does  not  affect 
the  modeling  of  the  heat  transfer  behavior. 

There  is  a  gas  turbine  parameter  relating  the  turbine  tip  speed 
to  the  speed  of  sound  in  the  working  fluids.  This  parameter  does  not 
appear  in  the  model  groups  as  the  differential  equations  are  written 
for  essentially  incompressible  flow.  The  effect  of  high  velocities 
is  to  produce  variations  in  the  fluid  properties  normal  to  the  flow 
direction.  For  the  model,  the  variations  are  small,  and  the  fluid 
properties  are  essentially  constant.  For  extension  to  high  velocity 
flow,  the  results  can  be  modified  using  the  temperature  dependent 
property  correlations  of  (14).  Therefore,  the  disk  peripheral  Mach 
number  is  not  a  significant  parameter  for  the  model  tests. 

A  survey  of  the  radial  flow  gas  turbine  field  reveals  the  fol¬ 
lowing  nondimensional  parameter  ranges  to  be  of  interest: 


nr,  f 

= 

105 

to 

106 

nr,m 

= 

105 

to 

5  x  106 

nr.m 

NP  F 

= 

CM 

O 

to 

1 

★ 

zo 

= 

0.1 

to 

0.3 

Npr 

= 

0.7 
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IV 


TEST  FACILITY 


A .  Description  of  Test  Facility 

The  present  model  is  a  rotating  disk  without  vanes  over  which 
air  flows  radially  inward.  The  disk  is  heated  relative  to  the  colder 
air  flow.  The  fluid  behavior  and  the  heat  transfer  coefficients  have 
been  determined  for  this  model.  An  overall  description  of  the  test 
facility  follows,  with  the  major  components  described  in  detail  in 
the  later  sections. 

A  schematic  of  the  test  facility  is  shown  in  Fig.  3,  and  the 
actual  test  facility  in  Fig.  4.  The  disk  is  driven  through  a  flat  belt 
pulley  drive  by  a  variable  speed  d-c  motor,  and  spins  about  a  verti¬ 
cal  axis.  The  disk  is  heated  by  the  radiant  energy  emitted  from  the 
tubular  quartz  infrared  heaters  (Fig.  8).  The  distribution  of  the 
radiant  energy  is  controlled  by  the  vane  reflectors  to  produce  a  uni¬ 
form  disk  test  surface  temperature.  The  radiant  energy  is  absorbed 
by  the  upper  surface  of  the  disk,  conducted  through  the  disk,  and 
convected  to  the  air  flowing  over  the  test  surface.  The  heat  flow 
through  the  laminated  aluminum  disk  is  measured  by  the  heat  meters 
bonded  in  the  transite  center  disk  (Fig.  5).  The  disk  test  surface 
temperature  is  measured  by  thermocouples  embedded  in  the  test  sur¬ 
face  (Fig.  6). 

The  air  stream  is  induced  over  the  test  surface  by  a  centrifugal 
blower.  Provisions  are  made  for  a  heated  veil- flow  over  the  test 
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surface.  The  screens  affixed  to  the  disk  periphery  (Fig.  6)  accelerate 
the  air  flows  to  the  disk  peripheral  velocity;  the  air  leaves  the  screens 
in  the  radial  direction  with  respect  to  the  disk.  The  air  mass  flow 
rate  is  measured  in  the  metering  section  upstream  of  the  blower  inlet. 

The  output  leads  of  the  heat  meters  and  thermocouples  mounted 
on  the  disk  lead  up  through  the  hollow  drive  shaft  to  the  rotating 
selector  switch  (Figs. 9,  10  and  11).  The  selector  switch  outputs  go 
to  the  garter  spring  pickup.  ^  The  pickup  transfers  the  output  volt¬ 
ages  from  the  rotating  shaft  to  a  stationary  p^+entiometer  for  measure¬ 
ment  of  the  heat  meter  and  thermocouple  voltages. 

The  air  stream  velocity  and  direction  are  measured  in  the  flow 
passage  between  the  rotating  disk  and  the  shroud.  A  three-hole  total 
pressure  and  yaw  probe  is  moved  across  the  flow  passage  by  a  trav¬ 
ersing  mechanism  (Fig.  13).  The  yaw  probe  angle  is  measured  by  a 
protractor  attached  to  the  traversing  mechanism.  The  fluid  static  pres¬ 
sure  is  obtained  from  static  pressure  taps  located  in  the  plexiglass 
shroud  (Fig  .  12) . 

The  heat  transfer  and  fluid  dynamics  parameters  are  determined 
from  these  measured  quantities.  The  disk  speed  is  measured  with  a 
Strobotac.  The  air  mass  flow  rate  is  measured  at  the  metering  sec¬ 
tion  by  the  A.S.M.E.  standard  orifice  installation.  The  fluid  velocity 

^Developed  on  the  suggestion  of  Mr.  W.  A.  Casler  of  Aerojet 
General  Corporation .  Patents  on  the  garter  spring  pickup  are 
pending . 
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and  flow  direction  are  computed  from  the  measured  angles  and  dynamic 
pressures.  The  air  properties  are  determined  from  temperature  and 
pressure  measurements  and  standard  tabulated  equation  of  state  infor¬ 
mation  (S) .  Heat  transfer  coefficients  are  evaluated  from  the  heat 
flux  measured  by  the  heat  meters  and  the  temperature  difference 
between  the  disk  test  surface  and  the  incoming  air. 

In  the  following  sections,  the  major  components  of  the  test 
facility  are  described  in  more  detail. 


B .  Rotating  Disk 

The  disk  design  was  governed  by  model  laws  and  heat  transfer 
considerations.  The  design  NR  M  operating  range  is  105  to  2  x  10  , 
a  range  representative  of  gas  turbine  operations.  A  disk  diameter 
of  16"  and  a  maximum  speed  of  8000  RPM  were  chosen  as  having  a 
low  operating  centrifugal  stress  and  of  convenient  size  for  construc¬ 


tion.  Aluminum  was  chosen  for  the  disk  material  as  the  stress 


level  is  the  lowest  and  the  factor  of  safety  the  highest  relative  to 
other  conventional  materials . 


The  disk  thickness  was  dictated  by  the  chosen  method  of 
determining  the  local  heat  transfer  coefficients.  The  local  heat  flux 
at  the  disk  test  surface  is  determined  from  the  heat  meter  surface 


heat  flux  and  the  test  surface  temperature  distribution  using  the  ana¬ 
lytical  solution  to  this  heat  conduction  problem  (Appendix  D) .  A 
uniform  test  surface  temperature  relative  to  the  temperature  difference 
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between  the  disk  end  air  was  desired.  The  final  disk  thickness  was 
1.682  inches. 

Preliminary  calibration  tests  selected  as  the  most  satisfactory 
heat  meters  those  manufactured  by  National  Instrument  Laboratories, 
Inc.,  Washington,  D.C.  The  heat  meter  consists  of  a  silver- 
tellerium  wafer  (0.4  inches  in  diameter  and  0.0  6  inches  thick)  with 
a  fine  mesh  copper  screen  bonded  to  each  side.  The  temperature 
difference  between  the  two  copper  and  silver-tellerium  junctions 
due  to  heat  flow  normal  to  the  wafer  produces  a  net  EMF  proportional 
to  the  heat  flux.  The  calibration  factors  of  these  meters  are  in  the 
range  600  -  800  (BTU/hr  ft^)/mv.  (A  detailed  description  of  the 
calibration  is  given  in  Appendix  E.) 

The  upper  surface  of  the  thick  aluminum  disk  was  instrumented 
with  five  heat  meters.  (See  Fig  5  for  a  detail  of  the  heat  meter 
placement  and  disk  construction.)  The  heat  meters  are  placed  at 
the  area-average  radius  of  each  of  five  equal  area  annular  sections. 
The  heat  meters  are  mounted  on  a  transite  disk,  a  material  of  about 
the  same  thermal  conductivity  as  the  heat  meters  .  This  insures  that 
the  thermal  resistance  normal  to  the  disk  will  be  uniform  over  the 
disk  surface.  To  protect  the  heat  meter  section,  a  thin  aluminum 
disk  was  bonded  over  the  transite  disk.  Copper-constantan  thermo¬ 
couples  were  mounted  in  this  disk  to  measure  the  heat  meter  tempera¬ 
tures.  The  disk  assembly  was  bonded  together  with  an  epoxy  resin 
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(Armstrong  A-12) .  The  heat  meters  . were  calibrated  after  bonding 
(Appendix  E) . 

The  disk  test  surface  was  instrumented  with  five  copper- 
constantan  thermocouples  placed  at  the  same  radii  as  the  heat 
meters.  The  test  surface  is  shown  in  Fig.  6.  The  thermocouples 
are  placed  in  grooves  milled  along  isotherms  (circumferential  lines) 
to  reduce  heat  conduction  along  the  leads .  The  thermocouple  end 
is  peened  into  the  groove  to  insure  good  thermal  contact  with  the 
disk.  The  leads  are  electrically  insulated  from  the  disk  by  a  coat 
of  Glyptal  cement. 

The  assembled  disk  specifications  are: 

Aluminum  Alloy  No.  2024  -  T4 

Yield  Stress  (0.2%  offset)  47,000  psi 

Maximum  Stress  (computed  at  8000  RPM)  9,000  psi 


Disk  Diameter  15.958 

Disk  Thickness:  Test  disk  1  •  682 

Heat  meter  section  0.0  63 

Epoxy  resin  0.006 

Protective  disk  0-065 

Total  thickness  1.816 


Thermal  Conductivity: 

Disk  (6)  72  BTU/(hr  ft2  F/ft) 

Heat  meter  section  0.37  BTU/(hr  ft2  F/ft) 
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C .  Inlet  Velocity  Control 


The  design  objective  for  the  inlet  velocity  control  was  to  pro¬ 
duce  a  flow  that  was  radial  with  respect  to  the  disk  at  the  disk 
periphery.  This  flow  condition  is  representative  of  radial  flow  gas 
turbine  operation.  Preliminary  tests  showed  that  a  simple  variable 
pitch  vane  system  installed  in  the  inlet  passage  could  not  turn  the 
flow  sufficiently  for  large  ratios  of  disk  tangential  velocity  to 
through  flow  velocity  (large  p) .  It  was  decided  that 

more  elaborate  methods  using  vanes  would  not  have  the  flexibility 
desired  for  the  test  facility. 

Preliminary  tests  also  showed  that  screens  wound  around  the 
disk  periphery  and  extending  over  the  air  inlet  passage  would  pro¬ 
duce  the  desired  flow  condition.  The  air  flow  would  be  accelerated 
by  the  screens  to  the  disk  tangential  velocity,  and  leave  the  screens 
in  the  radial  direction  relative  to  the  disk.  A  design  using  10  layers 
of  50  mesh  aluminum  screen  0.009  inches  in  diameter  was  chosen. 
Aluminum  was  used  because  of  its  low  mass  and  consequent  low 
centrifugal  forces.  A  fine  mesh  minimized  the  heat  transfer  from  the 
disk  along  the  screen  to  the  air  stream. 

A  continuous  4.5"  wide  strip  of  screen  was  wound  over  the  disk 
and  inner  lower  hoop.  The  two  outer  hoops  were  then  placed  over  the 
screen  and  bolted  to  the  disk  and  inner  hoop  respectively.  These 
hoops  are  made  of  4130  aircraft  steel  (98,000  psi  yield  stress) 
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to  withstand  the  centrifugal  loading  of  both  hoop  and  screen.  The 
resulting  air  inlet  passage  height  is  2.4  inches;  the  maximum  zQ* 
obtainable  is  0.3.  The  screens  installed  are  shown  in  Fig.  6.  The 
inlet  velocity  passage  is  shown  schematically  in  Fig.  7. 

The  disk  screen  assembly  was  balanced  on  a  commercial 
balancing  machine  to  4000  RPM .  The  unbalance  is  within  8  gm-in, 
corresponding  to  an  eccentricity  of  0.0005  inches.  The  runout  of 
the  lower  hoops  is  ±  0.050  inches.  The  natural  frequency  of  the 
disk  assembly  mounted  in  the  test  facility  is  in  the  range  2800  - 
3000  RPM. 

D .  Infrared  Heating  System 

The  design  requirement  for  the  disk  heating  system  was  that 
a  heat  flux  varying  with  radius  in  a  controlled  manner  could  be 
applied  to  the  disk.  Preliminary  tests  showed  that  the  infrared 
heating  method  would  satisfy  the  design  objective,  with  the  advan¬ 
tage  over  more  conventional  heating  methods  that  an  airtight  seal 
at  the  disk  periphery  would  not  be  necessary.  The  installed  heating 
system  is  shown  in  Fig.  3. 

Heating  is  accomplished  by  four  General  Electric  Tubular 
Quartz  heaters  producing  a  total  of  6400  watts.  The  reflectors  are 
gold-plated  copper  sheet;  gold-plating  is  needed  to  obtain  a  high 
reflectivity  in  the  infrared  spectrum.  A  seal  on  the  disk  drive  shaft 
minimizes  the  airflow  through  the  reflector  system. 
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The  radiant  energy  absorbed  at  the  upper  surface  of  the  disk  can 
be  controlled  in  three  ways.  The  heater  power  intensity  is  controlled 
by  varying  the  heater  voltage.  The  radiant  energy  distribution  imping¬ 
ing  on  the  disk  is  controlled  by  varying  the  vane  reflector  angle.  The 
energy  absorbed  by  the  disk  can  be  controlled  by  selectively  painting 
the  disk  with  a  high  absorbtivity  paint.  For  the  present  model,  the 
upper  surface  was  covered  with  a  uniform  coating  of  blackboard 
slating  with  an  absorbtivity  of  approximately  0.9.  The  radiant  energy 
distribution  absorbed  by  the  disk  was  thus  controlled  by  voltage  con¬ 
trol  and  vane  angle. 

E .  Selector  Switch  and  Garter  Spring  Pickup 

The  use  of  thermoelectric  elements  mounted  on  the  rotating  disk 
necessitated  the  development  of  a  selector  switch  and  pickup  assembly 
for  transmitting  the  electrical  outputs  to  stationary  recording  instru¬ 
ments.  The  unit  developed  for  this  purpose  is  depicted  in  Figs.  9  and 
10.  In  Fig.  11,  the  unit  is  shown  installed  on  the  disk  drive  shaft. 

The  30  leads  from  the  disk  instruments  (two  each  from  the  five  heat 
meters  and  two  each  from  ten  thermocouples)  are  brought  up  through 
the  hollow  disk  drive  shaft  to  the  copper  pins  of  the  connector  (Fig.  II)  . 
The  copper  leads  from  the  connector  sockets  are  soldered  to  the  copper 
selector  contacts  bonded  in  the  selector  switch  housing  (Fig.  9).  Four 
circular  segment  contacts  are  located  at  each  axial  position  to  accom¬ 
modate  four  leads  from  two  thermoelectric  elements.  The  selector  rod 
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is  moved  axially  to  obtain  outputs  from  any  one  of  the  eight  sets  of 
contacts  .  Copper  leads  (not  shown  in  Fig .  S)  from  the  copper  foil 
contacts  are  brought  through  the  hollow  selector  rod  to  the  rotating 
vee-groove  garter  spring  contacts  . 

The  outputs  are  transferred  from  the  rotating  to  the  stationary 
garter  spring  contacts  through  the  garter  springs.  These  springs 
are  made  from  beryllium  copper  wire  wound  to  form  a  helical  spring, 
with  the  ends  of  the  spring  soldered  together  to  form  a  circle.  The 
spring  rolls  between  the  rotating  and  stationary  vee-groove  contacts 
in  the  manner  of  a  planetary  gear  train.  Leads  from  the  four  sta¬ 
tionary  garter  spring  contacts  are  taken  to  a  potentiometer  for 
measurement  of  the  element  outputs. 

The  circuit  for  measurement  of  disk  thermocouple  outputs  is 
modified  to  account  for  the  constantan-copper  junctions  in  the  con¬ 
nector.  A  stationary  copper-constantan  compensation  junction  is 
introduced  in  the  circuit  to  oppose  the  connector  constantan-copper 
junctions.  The  compensation  junction  is  maintained  at  the  same 
temperature  as  the  connector  junctions  by  a  small  blower  that  con¬ 
tinuously  circulates  air  over  the  junctions.  The  ventilated  connec¬ 
tor,  compensation  junction,  and  blower  are  also  seen  in  Fig.  11. 

An  ice  junction  is  used  as  the  reference  junction  for  the  thermo¬ 
couple  output  measurement. 

The  primary  advantage  of  the  garter  spring  pickup  over 
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conventional  slip  ring  or  mercury  pool  units  is  the  lack  of  significant 
extraneous  EMF  production  at  high  rotative  speeds.  This  is  due  to 
the  absence  of  slipping,  with  consequent  rubbing  and  heating,  between 
the  spring  and  the  vee-grooves.  Test  showed  that  a  maximum  noise 
level  of  -  5  microvolts  (t  0.2  F)  was  produced  at  speeds  up  to  8000 
RPM.  Observations  of  the  spring  motion  with  a  Strobotac  verified 
the  lack  of  slippage,  and  showed  that  the  spring  remains  in  contact 
with  the  vee-grooves  at  all  speeds.  The  simplicity  of  construction 
and  operation  are  also  very  desirable  features. 

F  .  Flow  System 

An  airflow  system  was  designed,  constructed,  and  instrumented 
to  provide  metered  main  and  veil-flows  (Fig.  3).  A  suction  system 
with  a  13  stage  centrifugal  vacuum  cleaner  is  used  for  the  main  flow 
system  to  eliminate  the  large  plenum  chamber  that  would  be  needed 
in  a  pressure  system.  A  compressor  is  used  to  supply  the  veil-flow 
so  that  the  flowrate  may  be  varied  independently  of  the  main  flow. 

Both  flows  are  metered  by  calibrated  A.S.M.E.  standard  orifice 
installations . 

Preliminary  testing  revealed  that  the  main  airflow  pressure 
drop  through  the  test  section  became  prohibitively  high  at  high  disk 
speeds  (high  m/Nr  p) ,  with  a  resulting  decrease  in  flow  rate. 

This  was  attributed  to  the  fluid  kinetic  energy  dissipation  in  the 
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right  angle  turn  from  the  passage  between  the  disk  and  the  shroud  to 
the  main  flow  ducting.  At  the  higher  disk  speeds  the  tangential 
velocities  were  on  the  order  of  400  ft/sec  (Section  V) ,  and  pressure 
measurements  showed  that  all  of  the  kinetic  energy  was  dissipated. 

To  take  the  swirl  out  of  the  air  and  recover  some  kinetic 
energy  before  the  turn  into  the  ducting,  a  vaned  diffuser  was 
designed  and  installed.  The  diffuser  mounted  on  the  shroud  is. 
shown  in  Fig.  12.  With  the  diffuser,  the  flow  rate  at  the  maximum 
test  speed  was  doubled,  and  a  nearly  constant  flow  rate  was  obtain¬ 
able  over  the  test  range.  The  change  in  flow  pattern  affected  about 
10  per  cent  of  the  disk  heat  transfer  area.  However,  the  effect  on 
the  mean  heat  transfer  results  is  insignificant  (Section  VII). 

G .  Fluid  Dynamics  Instrumentation 

The  fluid  dynamics  was  determined  from  measurements  of  the 
fluid  velocity,  flow  angle,  and  wall  static  pressure.  The  wall 
static  pressure  was  measured  through  wall  taps  located  on  the 
shroud  (Fig.  12).  As  the  pressure  tap  diameter  was  fairly  large 
(0.128  inches),  inserts  were  made  to  increase  the  length  of  the 
pressure  tap  hole  to  2 . 5  diameters  in  accordance  with  the  A.S.M.E. 
recommendations  (1)  .  The  inserts  were  removed  when  the  probe  was 
entered  through  the  tap  for  measurement  of  the  fluid  velocity. 

The  fluid  velocity  and  flow  angle  were  measured  by  a  two- 
dimensional  directional  probe  manufactured  by  the  United  Sens  or 
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and  Control  Corporation,  Glastonbury,  Conn.  (Fig.  13).  The  probe 
tip  consists  of  three  parallel  open-ended  tubes  0.0  30  inches  in 
diameter.  The  two  side  tubes  are  chamfered  30p,  and  all  tubes 
are  bent  90°  so  as  to  measure  pressures  perpendicular  to  the 
probe  axis,  as  shown  in  the  following  figure.  The  center  tube  is 
used  to  measure  the  fluid  total  pressure.  The  dynamic  pressure  is 
measured  as  the  difference  between  the  total  pressure  and  the 
wall  static  pressure.  The  fluid  resultant  velocity  was  computed 
from  this  dynamic  pressure. 

The  pressures  in  the  two  side  tubes  were  balanced  by  probe 
rotation  to  determine  the  flow  angle.  A  protractor  affixed  to  the 
traversing  mechanism  allowed  measurement  of  the  amount  of  rota¬ 
tion.  The  fluid  velocity  components  were  computed  from  the 
resultant  velocity  and  flow  angle. 


Top  View  I - ^ 


Air  Flow 


Side  View 
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Air  Flow 


Probe  Axis 
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FIGURE  3 


SCHEMATIC  OF  ROTATING  DISK  TEST  FACILITY 


FIGURE  4 


ROTATING  DISK  TEST  FACILITY 
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SCHEMATIC  OF  DISK  CROSS  SECTION  AND 
INSTRUMENT  LOCATIONS 

F IGURE  6 

DISK  TEST  SURFACE  AND  ROTATING  SCREENS 

FIGURE  7 

SCHEMATIC  OF  TEST  FACILITY  INLET  PASSAGE 

FIGURE  8 

INFRARED  HEATING  SYSTEM 
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FIGURE  9 


SCHEMATIC  OF  THE  SELECTOR  SWITCH  AND 
GARTER  SPRING  PICKUP  ASSEMBLY 

1 .  Switch  and  Pickup  Support 

2.  Selector  Switch  Housing 

3.  Selector  Switch  Contacts 

4.  Copper  Foil  Contacts 

5 .  Selector  Rod 

6.  Garter  Spring  Housing 

7 .  Garter  Spring 

8.  Garter  Spring  Contacts 

9 .  Ball  Bearing 

10.  Guide  Rod  for  Sliding  the  Garter  Spring  Housing 


FIGURE  10 

CROSS  SECTION  THROUGH  SELECTOR  SWITCH 
AND  GARTER  SPRING  PICKUP 

FIGURE  1 1 

SELECTOR  SWITCH  AND  GARTER  SPRING  PICKUP 
MOUNTED  ON  DISK  DRIVE  SHAFT 
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SECTION  THROUGH  GARTER 
SPRING  PICKUP 


FIGURE  12 


SHROUD  AND  VANED  DIFFUSER 


FIGURE  13 

DIRECTIONAL  PROBE  AND  TRAVERSING 
MECHANISM 
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SEALING  SEC  TION 
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PROBE  T  I  P 


V  EXPERIMENTAL  RESULTS 


A .  Radial  Inflow  Flow  Results 

The  objectives  of  the  flow  investigation  are  to: 

1.  Evaluate  the  performance  of  the  rotating  screens  and 
sealing  section . 

2.  Experimentally  determine  the  fluid  dynamics  behavior. 

3.  Compare  the  experimentally  determined  fluid  dynamics 
behavior  with  the  behavior  predicted  by  the  analysis 
(Section  VI) . 

The  component  performance  was  evaluated  during  the  radial 
inflow  tests.  Therefore,  the  experimental  flow  results  will  be  pre¬ 
sented  so  as  to  answer  the  objectives  concurrently. 

The  flow  behavior  was  studied  at  a  single  disk- to- shroud 
clearance  ratio  of  z0*  =  0.113  over  a  machine  Reynolds  number  range 
Nd  M  =  3.0xl04to7.5xl05  and  a  flow  Reynolds  number  range 
Nr  F  =  3.4  x  105  to  5.8  x  105.  The  test  NRjM/NRjF  range  was  0.05 
to  1.6.  Four  types  of  flow  results  were  obtained: 

1.  Tangential  and  radial  component  velocity  profiles  and  flow 
angle  profiles  were  determined  at  five  radial  positions  tor 
five  Nj^m/Nr^p  ratios. 

2.  Total  velocity  and  flow  angle  profiles  were  determined  at 
four  circumferential  positions  for  two  radial  positions 
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ratios . 


and  four  Nr  ,m/Nr  f 

3.  Total  velocity  and  flow  angle  profiles  were  determined  at 
two  radial  positions  for  four  to  six  values  of  Nr  p  at  two 
Nr  m/Nr  f  ratios  •  This  was  done  to  experimentally  verify 
that  the  nondimensional  parameters  correlated  the  fluid 
dynamics  behavior. 

4.  Wall  static  pressure  variations  with  radial  position  were 
measured  at  four  circumferential  positions  for  five 

nr,m/nr,f  ratios- 

The  fluid  velocity  and  flow  angle  was  measured  with  the  two- 
dimensional  directional  probe.  The  tangential  and  radial  components 
were  computed  from  the  flow  angle  and  the  resultant  velocity.  The 
wall  static  pressure  was  measured  with  the  static  pressure  taps 
located  in  the  shroud. 

The  radial  and  tangential  velocity  profiles  and  the  flow  angle 
profiles  for  five  radial  positions  are  presented  in  Figs.  14  to  18  for 
five  Nr  m/Nr  p  ratios.  The  corresponding  profiles  predicted  from 
the  one-dimensional  flow  analysis  (Analytical  Investigation,  Sec¬ 
tion  VI,  Eq.  [22] s  [  231  and  [24]  )  are  also  indicated  on  the 
figures.  The  profiles  are  predicted  to  be  uniform  across  the  flow 
passage . 

The  measured  profiles  at  r*  =  0.893  show  that  the  entering 
velocity  profile  is  not  uniform  as  desired.  The  shape  of  the  profile 
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results  from  flow  leakage  around  the  edge  of  the  rotating  screen  in  the 
sealing  section  (Fig.  7).  The  leakage  flow  enters  the  main  flow  passage 
along  the  screens  ,  creating  a  stalled  region  around  the  circumference 
of  the  plexiglass  shroud.  As  the  leakage  flow  has  not  been  accelerated 
to  the  disk  tangential  velocity  by  the  screens,  a  backflow  in  both  the 
radial  and  tangential  directions  is  created  in  the  stalled  region. 
Schematically,  the  flow  picture  near  the  test  section  inlet  is: 


The  extent  of  the  stalled  region  is  greatest  at  low  NRjM/NR(F 
(0.056,0.28).  At  high  NR  M/NR  p  (1.20,1.56),  the  effect  of  the 
leakage  flow  is  considerably  less;  the  stalled  region  does  not  extend 
as  far  radially  inward  and  the  entering  profiles  are  more  uniform. 

The  inlet  profile  appears  to  affect  the  flow  pattern  only  near 
the  outer  edge  of  the  disk.  This  is  especially  evident  in  Fig.  16, 
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where  the  r*  =  0.893  profile  is  completely  unlike  the  shape  of  the 
other  profiles,  while  the  r*  =  0 . 768  profile  is  very  similar  to  the 
remaining  profiles.  It  appears  that  the  effects  of  the  disk,  shroud, 
and  flow  area  decrease  outweigh  the  effect  of  the  entering  profile. 

The  effect  of  the  entering  profile  is  further  lessened  at  higher 

nr,m/nr,f° 

A  striking  feature  of  the  flow  pattern  is  that  at  high  Np  p 

the  air  is  flowing  radially  inward  near  the  disk  and  shroud,  yet 
outward  near  the  passage  center.  As  seen  in  Figs.  17  and  18  for 
Nr  M/Nr  f  1-20  and  1.56  respectively,  the  radial  outflow  velocity 
is  about  equal  to  the  radial  inflow  velocity  near  the  disk,  and  one-fifth 
that  near  the  shroud.  This  anomalous  behavior  was  not  predicted  by 
the  analysis.  A  three-dimensional  representation  of  a  profile  illus¬ 
trating  this  radial  outflow  is  presented  in  Fig.  19. 

The  existence  of  backflow  in  the  passage  center  can  be 
rationalized.  As  the  air  flows  radially  inward,  the  angular  momentum 
of  the  air  is  essentially  conserved  (as  shown  in  Fig.  21  and  discussed 
later).  Consequently,  the  air  tangential  velocity  increases  in  the 
radially  inward  direction.  The  disk  tangential  velocity  decreases  in 
the  radially  inward  direction,  and  is  less  than  the  fluid  tangential 
velocity.  As  a  result,  the  wall  shearing  stresses  on  both  the  disk 
and  the  shroud  act  to  decrease  the  air  tangential  velocity  and  create 
boundary  layers  of  low  tangential  velocity  fluid. 
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The  centrifugal  force  field  acting  on  the  air,  proportional  to  the 
square  of  the  tangential  velocity,  is  consequently  greatest  in  the  flow 
passage  center.  At  high  j^/Np  f<  this  force  is  evidently  high 
enough  to  force  the  fluid  itself  to  flow  radially  outward.  The  lower 
centrifugal  force  field  in  the  regions  near  the  disk  and  shroud  allow 
the  fluid  to  flow  radially  inward.  As  the  disk  is  rotating,  the  centrifu¬ 
gal  force  is  higher  near  the  disk  than  the  shroud.  Consequently,  the 
radial  inflow  in  the  region  near  the  shroud  is  greater  than  that  near 
the  disk . 

In  Fig.  20,  a  cross  section  of  the  flow  pattern  is  shown  for 
Np  m/NR  F  =  1-20.  This  flow  pattern  is  somewhat  similar  to  the 
flow  pattern  created  by  water  draining  from  a  tank.  In  this  case, 
the  angular  momentum  of  the  water  is  conserved  as  the  water  flows 
radially  inward.  At  some  radius,  the  resulting  centrifugal  force  is 
high  enough  to  prevent  further  radial  inflow,  resulting  in  a  hollow 
vortex  core . 

The  flow  rates  obtained  from  the  integrated  radial  velocity  pro¬ 
files  were  compared  to  the  flow  rate  measured  at  the  flow  metering 
section.  The  difference  was  within  the  accuracy  of  determining  the 
radial  velocity  profiles.  The  uncertainty  in  the  profiles  is  very  high 
at  large  Nr,m/Nr  p  (the  velocity  uncertainty  ranges  are  1  0.5  per 
cent  at  Np(  M/NR  p  =  0  . 056,  t  6.5  per  cent  at  0 . 61,  and  ±  20  per 
cent  at  1.5  7).  The  flow  angle  is  near  90°,  and  the  radial  velocity 
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is  computed  as  the  cosine  of  the  flow  angle  times  the  measured 
resultant  velocity.  Consequently,  small  errors  in  flow  angle  produce 
large  errors  in  the  radial  velocity  component.  Nevertheless,  the 
agreement  between  integrated  and  metered  flows,  and  the  self- 
consistency  of  the  velocity  profiles,  strengthens  the  experimentally 
determined  picture  of  the  flow  pattern. 

The  tangential  velocity  profiles  are  consistently  lower  than  the 
profiles  predicted  by  the  analysis.  This  is  a  result  of  flow  bypassing 
the  rotating  screens.  As  a  measure  of  the  performance  of  the  screens, 
the  angular  momentum  flux  for  the  four  highest  NR  ^/NR  p  ratios  was 
evaluated.  The  angular  momentum  flux  is  related  to  the  shearing 
stresses  on  the  disk  and  shroud  by  the  following  equation  (Analytical 
Investigation,  Section  VI): 


( t 


6  disk 


6  shroud^ 


dr* 


9c 


zo  ro  r 


u  v  d  zi 


[9] 


where 


angular  momentum  flux  = 


do] 


The  angular  momentum  flux  is  presented  as  a  function  of  radius 
for  the  four  Nr(m/nR,F  ratios  in  Fig.  21a.  The  flux  is  constant 
within  the  experimental  uncertainty,  showing  that  the  disk  and 
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shroud  shear  stresses  are  small  relative  to  the  momentum  fluxes. 

The  constant  value  of  momentum  flux  provides  a  check  on  the  instru¬ 
mentation.  Any  significant  increase  in  momentum  flux  in  the  radially 
inward  direction  would  violate  the  laws  of  motion,  and  therefore  be  a 
result  of  inaccuracies  in  measuring  velocity. 

The  ratio  of  the  actual  angular  momentum  flux  to  that  predicted 
is  plotted  as  a  function  of  NR<M/NR  p  in  Fig.  21b.  The  values  of 
0.70  -  0.75  indicate  that  25  -  30%  of  the  main  flow  bypasses  the 
rotating  screens,  and  is  not  accelerated  to  the  disk  tangential 
velocity.  The  apparent  increase  in  the  ratio  could  be  due  to  the 
experimental  inaccuracy  of  determining  the  radial  velocity  component. 
The  uncertainty  is  highest  at  high  NR,jyj/NR  p,  and  the  integrated 
profiles  tended  to  yield  higher  flow  rates  than  those  measured  at  the 
orifice  metering  section.  However,  an  increase  in  the  momentum 
flux  ratio  with  increased  NR^M/NR,p  is  supported  by  comparing  the 
measured  pressure  drop  across  the  screens  to  the  pressure  drop  com¬ 
puted  for  no  bypass  flow;  these  results  indicate  that  less  bypass 
flow  occurs  at  higher  disk  speeds  (higher  NR  m/N  ). 

The  variation  of  velocity  and  flow  angle  profile  with  circum¬ 
ferential  position  for  four  NR  jyj/Np  p  ratios  is  shown  in  Figs.  22  to 
25.  In  Fig.  2  6,  the  tangential  velocity  variation  with  circumferential 
position  is  shown  for  two  NR^/NR  p.  Axial  symmetry  is  seen  to  be 
best  at  the  higher  NR  jyj/NR  p.  The  lack  of  axial  symmetry  probably 
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results  from  the  eccentricity  of  the  sealing  section  with  respect  to 
the  disk  axis.  The  bypass  flow  area,  and  therefore  the  bypass  flow, 
would  then  vary  around  the  shroud  periphery.  This  would  result  in 
an  asymmetric  main  flow.  It  is  doubtful  that  the  asymmetry  was 
caused  by  variations  in  the  shroud-to-disk  clearance,  as  the  varia¬ 
tion  in  clearance  around  the  periphery  was  less  than  ±  0.  7  percent. 

in  conclusion,  at  NR  j^/Nj^  p  of  0.61  and  greater,  the  flow 
is  essentially  axisymmetric .  This  is  consistent  with  the  previous 
velocity  profile  results,  where  the  effect  of  the  inlet  is  less  at 
higher  NRjM/NR  p  . 

Nondime nsional  velocity  and  flow  angle  profiles  were  compared 

to  help  verify  that  the  nondimensional  parameters  chosen  govern  the 

fluid  dynamics.  In  Figs.  2  7  and  28,  comparisons  are  made  with  Nn  r 

k,  r 

as  a  parameter  for  various  Nr>m/Nr^f  ratios.  The  nondimensional 
profiles  are  all  similar  in  shape.  The  discrepancy  between  the  pro¬ 
files  exceeds  the  estimated  uncertainty,  and  does  not  vary  systemati¬ 
cally  with  NR,p.  It  is  possible  that  additional  errors  in  velocity 
measurement  are  introduced  when  NR  p  is  a  variable.  (See  Appendix 
B  for  a  discussion  of  velocity  measuring  errors.)  The  similarity  in 
profile  shape  supports  the  choice  of  nondimensional  parameters,  but 
a  definite  verification  awaits  further  tests  . 

The  flow  angle  profiles  of  Figs  .  2  7  and  28  are  in  good  agree¬ 
ment  with  the  exception  of  the  results  for  r*  =  0.768  and  Nr<m/Nr  p  = 
0.30.  The  discrepancies  exist,  and  are  probably  caused  by  increased 
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flow  leakage  past  the  rotating  screens  at  low  flows  (NR  p) .  Again, 
the  inlet  effect  is  mitigated  at  the  inner  radii  by  the  flow  area  decrease. 
At  the  higher  NR  m/Nr  r,  the  effects  of  the  disk,  shroud,  and  area 
change  predominate.  This  agreement  between  the  flow  angle  profiles 
supports  the  choice  of  nondimen sional  parameters. 

The  variation  of  the  wall  static  pressure  with  radius  is  pre¬ 
sented  in  Fig.  29.  Data  was  obtained  at  all  r*  for  only  one  circum¬ 
ferential  position,  and  supplemented  by  some  data  at  other  circum¬ 
ferential  positions.  The  axial  symmetry  as  indicated  by  the  static 
pressures  is  not  as  good  as  that  indicated  by  the  velocity  profile 
results.  This  is  partially  a  result  of  small  uncertainties  in  the  pres¬ 
sure  tap  location.  The  tap  location  uncertainty  due  to  shroud  eccen¬ 
tricity  with  the  disk  axis  is  about  -  0.080  inches,  or  it  0.01  rQ. 

At  small  radii,  the  pressure  variation  with  radius  follows  the 
trend  of  the  predicted  variation  (Analytical  Investigation,  Section  VI, 
Eq.  [  25]  ).  The  discrepancies  at  the  disk  inlet  are  due  to  the  stall 
created  by  the  leakage  flow.  The  inlet  effect  is  again  felt  only  at 
the  outer  radii,  and  is  less  at  higher  NR  j^/Nr  R . 

Some  limited  nondimensional  static  pressure  variations  for  two 
Nr  m/Nr  p  are  presented  in  Figs.  29c  and  29d.  The  correlation  is 
quite  good.  This  helps  support  the  choice  of  p  and  NR  j^/Nr,  p 
as  the  nondimensional  parameters  governing  the  fluid  dynamics. 
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The  following  conclusions  can  be  drawn  from  the  radial  inflow 


fluid  dynamics  investigation: 

1.  The  performance  of  the  rotating  screen  and  seal  is  less  than 

desired.  The  rotating  screens  accelerate  only  70  -  75  per 
cent  of  the  main  flow  to  the  disk  tangential  velocity.  The 
remaining  25  -  30  per  cent  of  the  flow  bypasses  the  rotating 
screens.  To  accurately  establish  the  flow  picture,  this 
defect  must  be  remedied. 

2.  The  flow  that  bypasses  the  rotating  screens  affects  the  flow 

behavior  only  near  the  edge  of  the  disk.  The  inlet  effect 
diminishes  as  NR  M/NR  F  increases. 

3.  The  fluid  dynamics  behavior  is  described  in  Figs.  14-18, 

and  differs  considerably  from  that  predicted  by  the  analysis 
(Section  VI).  The  reasons  for  the  discrepancies  between 
experiment  and  analysis  are  discussed  in  Section  VII. 

4.  The  radial  inflow  velocities  are  greatest  in  the  regions  near¬ 

est  the  disk  and  the  shroud.  At  NR  M/NR  p  0 . 61  and 
higher,  the  flow  is  radially  outward  in  the  passage  center. 

5.  Within  the  experimental  uncertainties,  the  fluid  angular 

momentum  is  conserved  as  the  fluid  flows  radially  inward. 

6.  The  predicted  nondimensional  parameters,  Np  p  and 

Nr,m/Nr  p,  appear  to  govern  the  fluid  dynamics.  A  defi¬ 
nite  verification  awaits  additional  tests. 
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FIGURE  14 


FLUID  VELOCITY  COMPONENT  AND  FLOW 
ANGLE  PROFILES 

NR  F  =  5.6xl05 

NR,M  =  3.2  x  104 

NR,m/nr,f  =0-056 

(a)  Tangential  Velocity  Profiles 

(b)  Radial  Velocity  Profiles 

(c)  Flow  Angle  Profiles 


DATA 


O  r*  -  0.893 
□  r*  =  0.768 
A  r*  =  0.642 
O  r*  =  0.517 
V  r*  =  0.  391 


•  Disk  tangential  velocity  at  r*  =  0.893 
▼  Disk  tangential  velocity  at  r*  =  0.391 


ANALYSIS  (SECTION  VI) 

Tangential  velocity,  v,  predicted  by  Eq.  [  23  ] 
Radial  velocity,  u,  predicted  by  Eq.  t  22  ] 
Flow  angle  predicted  by  Eq .  [  24  1 

- O,  a,  a,o,v 
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FIGURE  15 


FLUID  VELOCITY  COMPONENT  AND  FLOW 
ANGLE  PROFILES 


nr,f 

=  5.6  x  105 

nr,  m 

=  1.59  x  10' 

nr,m//nr,  f 

=  0.28 

(a) 

Tangential  Velocity  Profiles 

(b) 

Radial  Velocity  Profiles 

(c) 

Flow  Angle  Profiles 

DATA 

G 

r* 

=  0.893 

□ 

r* 

=  0.768 

A 

r* 

=  0.642 

<3> 

r* 

=  0.517 

V 

r* 

=  0.391 

• 

Disk  tangential  velocity  at  r*  =  0.893 

▼ 

Disk  tangential  velocity  at  r*  =  0 . 391 

ANALYSIS  (SECTION  VI) 

Tangential  velocity,  v,  predicted  by  Eq.[  23  ] 
Radial  velocity,  u,  predicted  by  Eq,  [  22] 
Flow  angle  predicted  by  Eq.  [  24  ] 

- 0,0,  A,0,  V 
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FIGURE  16 

FLUID  VELOCITY  COMPONENT  AND  FLOW 
ANGLE  PROFILES 


Np^  p  =  5.6x10^ 

Mr.M  =  3.4  Xio4 


N 


R,m/NR,F 


=  0.61 


(a)  Tangential  Velocity  Profiles 

(b)  Radial  Velocity  Profiles 

(c)  Flow  Angle  Profiles 
DATA 

©  r*  =  0.893 

□  r*  =0.768 

A  r*  =0.642 

0  r*  =  0.517 

V  r*  =  0.391 

•  Disk  tangential  velocity  at  r*  =  0.893 
▼  Disk  tangential  velocity  at  r*  =  0.391 

ANALYSIS  (SECTION  VI) 

Tangential  velocity,  v,  predicted  by  Eq.  [  23  ] 
Radial  velocity,  u,  predicted  by  Eq.  [  22  ] 
Flow  angle  predicted  by  Eq.  [  24  ] 

- O,  □,  A  ,0,V 
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FLUID  VELOCITY  COMPONENT  AND  FLOW 
ANGLE  PROFILES 


Nr  f  =  5.2  x  105 

Nr,m  =  6.2  x  10^ 

NR,m/NR,F  =  1,20 


(a) 

Tangential  Velocity  Profile 

(b) 

Radial  Velocity  Profile 

(c) 

Flow  Angle  Profile 

DATA 

0  r* 

=  0.893 

0  r* 

=  0.768 

A  r* 

=  0.642 

0  r* 

=  0.517 

V  r* 

=  0.391 

•  Disk  tangential  velocity  at  r*  =  0.893 
T  Disk  tangential  velocity  at  r*  =  0.391 

ANALYSIS  (SECTION  VI) 


Tangential  velocity,  v,  predicted  by  Eq.  [23] 
Radial  velocity,  u,  predicted  by  Eq.  [  22  ] 
Flow  angle  predicted  by  Eq.  [24  ] 

- 0,0,  A  ,  O,  V 
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FLUID  VELOCITY  COMPONENT  AND  FLOW  ANGLE 
PROFILES 


NR,F 

=  4.8  x  105 

nr,m 

=  7.5  x  10: 

NR,m/NRjF 

=  1.56 

(a)  Tangential  Velocity  Profiles 

(b)  Radial  Velocity  Profiles 

(c)  Flow  Angle  Profiles 


DATA 


0 

r* 

=  0.893 

0 

r* 

=  0.768 

A 

r* 

=  0.642 

<•> 

r* 

=  0.517 

V 

r* 

=  0.391 

•  Disk  tangential  velocity  at  r*  =  0.893 
▼  Disk  tangential  velocity  at  r*  =  0.391 

ANALYSIS  (SECTION  VI) 

Tangential  velocity,  v,  predicted  by  Eq.  [23] 
Radial  velocity,  u,  predicted  by  Eq.  [  22  ] 
Flow  angle  predicted  by  Eq.  [  24  ] 
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FIGURE  19 

THREE-DIMENSIONAL  VELOCITY  PROFILE 
REPRESENTATION 


FIGURE  20 

FLOW  PATTERN  ILLUSTRATING  BACKFLOW 

/ 

IN  THE  CENTER  OF  THE  FLOW  PASSAGE 
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FIGURE  21 

FLUID  ANGULAR  MOMENTUM  FLUX 

(a)  Angular  Momentum  Flux  as  a  Function 
of  Radius 

(b)  Ratio  of  Actual  Momentum  Flux  to  that 
Predicted  as  a  Function  of  NR  M/NR  p 
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FIGURE  22 


RESULTANT  VELOCITY  AND  FLOW  ANGLE  PROFILE 
VARIATION  WITH  CIRCUMFERENTIAL  POSITION 


NR,F 

=  5. 6  x  10 

nr,m 

=  3.2  x  10 

nr,m/nr,f 

=  0.056 

(a)  Resultant  Velocity,  r*  =  0.768 

(b)  Flow  Angle,  r*  =  0.768 

(c)  Resultant  Velocity,  r*  =  0.517 

(d)  Flow  Angle,  r*  =  0.517 
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FIGURE  23 

RESULTANT  VELOCITY  AND  FLOW  ANGLE  PROFILE 
VARIATION  WITH  CIRCUMFERENTIAL  POSITION 


NR,F 

=  5. 6  x  105 

nr,m 

=  1.59  x  10! 

nr,m/nr,f 

=  0.28 

(a)  Resultant  Velocity,  r*  =  0.768 

(b)  Flow  Angle,  r*  =  0.768 

(c)  Resultant  Velocity,  r*  =  0.517 

(d)  Flow  Angle,  r*  =  0.517 
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FIGURE  24 


RESULTANT  VELOCITY  AND  FLOW  ANGLE  PROFILE 
VARIATION  WITH  CIRCUMFERENTIAL  POSITION 

NR/F  =  5.6  x  105 

Nr,M  -  3.4xl05 

NR,m/Nr,F  =  0  •  61 

(a)  Resultant  Velocity,  r*  =  0.768 

(b)  Flow  Angle,  r*  =  0.768 

(c)  Resultant  Velocity,  r*  =  0.517 

(d)  Flow  Angle,  r*  =  0.517 
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FIGURE  25 


RESULTANT  VELOCITY  AND  FLOW  ANGLE  PROFILE 
VARIATION  WITH  CIRCUMFERENTIAL  POSITION 

NR  p  =  5.2  x  105 

NR,M  =  6.2  x  105 

NR,M'/NR,F  =  1-20 

(a)  Resultant  Velocity,  r*  =  0.768 

(b)  Flow  Angle,  r*  =  0.768 

(c)  Resultant  Velocity,  r*  =  0.517 

(d)  Flow  Angle,  r*  =  0.517 
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FIGURE  2  6 


TANGENTIAL  VELOCITY  VARIATION  WITH 
CIRCUMFERENTIAL  POSITION 


(a)  nr.m/nr,f 

=  0.056 

(b)  nR-M/nRiF 

=  1.20 
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NONDIMENSIONAL  VELOCITY  AND  FLOW  ANGLE 
PROFILE  VARIATION  WITH  NR  p 


nr,m//nr,f 


=  0.30 


(a)  Nondimensional  Velocity,  r*  =  0.768 

(b)  Flow  Angle,  r*  =  0.768 

(c)  Nondimensional  Velocity,  r*  =  0.517 

(d)  Flow  Angle,  r*  =  0.517 
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FIGURE  28 

NONDIMENSIONAL  VELOCITY  AND  FLOW 
ANGLE  PROFILE  VARIATION  WITH  NR  p 


nr,m/nr,f 


=  1.22 


(a)  Nondimensional  Velocity,  r*  ■=  0.7  68 

(b)  Flow  Angle,  r*  =  0.768 

(c)  Nondimensional  Velocity,  r*  =  0.517 

(d)  Flow  Angle,  r*  =  0.517 
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FLOW  ANGLE  (degrees!  FLOW  ANGLE  (degrees) 


FIGURE  2  9 

WALL  STATIC  PRESSURE  VARIATION  WITH 
RADIAL  POSITION 

(a)  Static  Pressure,  ^/NR  p  -  0.045, 

0.28,  and  0 . 59 

(b)  Static  Pressure,  NR  M/NR  p  =  1.20 

and  1 . 58 

(c)  Nondimensional  Static  Pressure, 

nr,m/nr,f  =  0  -29 

(d)  Nondimensional  Static  Pressure, 

NR,m/Nr,  F  =  K19 
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B.  Radial  Inflow  Heat  Transfer  Results 


The  objectives  of  the  heat  transfer  study  are  to: 

1 .  Evaluate  the  technique  for  determining  the  local  and  average 

heat  transfer  coefficients. 

2.  Experimentally  determine  the  heat  transfer  behavior. 

3 .  Compare  the  experimentally  determined  heat  transfer 

behavior  with  the  predicted  behavior. 

The  experimental  technique  was  evaluated  during  radial  inflow 
heat  transfer  tests.  Therefore,  the  results  will  be  presented  so  as  to 
answer  the  objectives  concurrently. 

The  heat  transfer  behavior  was  studied  at  a  single  shroud- to-disk 
clearance  ratio  of  zQ*  =  0.113  over  a  flow  Reynolds  number  range  NR  p  = 
1.5  x  105  to  6.6  x  10s,  and  over  a  Reynolds  number  ratio  range 
NrjM/Nr  p  =  0.1  to  1.7.  Average  heat  transfer  coefficients  were 
calculated.  Local  heat  transfer  coefficients  could  not  be  satisfactorily 
determined . 

The  average  heat  transfer  behavior  is  represented  by  the  average 
Stanton  number,  Ngt  .  The  following  procedure  was  used  to  compute  N 

St 

The  average  heat  flux  through  the  disk  was  taken  as  the  area  average 
of  the  heat  fluxes  measured  by  the  heat  meters.  The  temperature  poten¬ 
tial  for  heat  transfer  was  taken  as  the  difference  between  the  area 
average  disk  surface  temperature  and  the  air  inlet  temperature.  The 
mass  velocity  was  determined  from  the  air  mass  flow  rate  and  the  test 
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section  inlet  area,  2  ti  rQ  zQ.  The  average  Stanton  number  is  thus 
defined  as: 
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St 


G  CL 


q"  /  (t 


disk 
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air  in' 
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2  it  rG  zQ 
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The  local  Stanton  numbers  Ngt  ^ ,  were  determined  using  the 
calculated  local  heat  flux  distribution  at  the  disk  test  surface  which 
is  determined  from  the  measured  heat  meter  heat  flux  and  the  test 
surface  temperature  distribution  using  Eq.  [D  35  ]  ,  Appendix  D: 

Jo  (  “m  rG  ) 
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where: 

q_"  (~  )  is  the  test  surface  heat  flux  distribution 

o  rQ 

tQ  ( n  )  is  the  test  surface  temperature  distribution 

qb"  (n  )  is  the  heat  meter  heat  flux  distribution 

The  Burroughs  220  Digital  Computer  at  Stanford  was  employed  to  cal¬ 
culate  q0"  )  from  Eq.  [  12 ]  . 

In  the  evaluation  of  the  local  Stanton  numbers,  the  definitions 
of  temperature  difference  and  mass  velocity  are  the  same  as  those 
for  the  average  Stanton  numbers.  The  local  Stanton  number  is  thus 
given  by: 
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The  average  heat  transfer  data  and  results  are  tabulated  in 
Table  I.  The  results  are  grouped  according  to  the  nominal  test  value 
of  nR,m/nr  p •  Runs  1  "  11  and  IS  -  22  are  data  taken  during  a 
preliminary  evaluation  of  the  test  facility.  Following  this  evalua¬ 
tion,  a  labyrinth  seal  was  added  to  the  inlet  section  to  reduce  the 
amount  of  flow  bypassing  the  rotating  screens  and  the  vaned  dif¬ 
fuser  was  installed.  The  reason  for  including  this  earlier  data  will 
be  discussed  later. 

The  heat  transfer  results  are  presented  in  Fig.  30.  The  co- 
0  2 

ordinates,  NstNR  p  *  and  NR  m/Nr  F'  are  those  predicted  by 

the  analysis.  The  data  is  not  correlated  by  these  parameters,  as 

0.2  . 

can  be  seen  by  the  spread  in  NgtNR  p  at  any  given  Nr^/Nr  p. 
The  data  appears  to  follow  the  trend  of  the  analysis,  showing 
increased  heat  transfer  as  NR fR/[/NR  p  increases.  However,  the 
analysis  appears  50  to  100  per  cent  low. 

Correlation  of  the  data  is  obtained  using  a  different  exponent 
of  the  flow  Reynolds  number.  This  exponent  is  determined  from 
cross  plots  of  Ngt  as  a  function  of  Nr  p  for  a  given  nominal  value 
of  Nr  m/Nr  F •  The  data  are  presented  in  this  manner  in  Fig.  31a 
and  Fig.  31b.  The  data  of  runs  1-4  and  IS  -  22  were  not  included 
in  the  curve  fitting.  It  was  felt  that  the  difference  in  inlet  geometry 
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between  the  preliminary  runs  and  the  remainder  of  the  data  produced 
a  significantly  different  flow  pattern  and  precluded  a  comparison  of 
the  data.  However,  it  is  very  interesting  to  note  that  the  high 
Nr  m/Nr  f  data  of  runs  5-11  are  compatible  with  the  later  data. 
This  reinforces  the  flow  investigation  conclusions  that  the  effect 
of  the  inlet  decreases  as  Nr  m/Nr  p  increases. 

Least  squares  straight  lines  were  passed  through  the  data  of 
Fig.  31.  The  exponent  of  Nr  p  (equal  to  the  negative  of  the  slope) 
is  tabulated  below: 


m/NR,F 

M  (-slope) 

0.15 

0.598 

0.31 

0.592 

0.50 

0.591 

0.66 

0.536 

0.S9 

0.546 

1.27 

0.525 

1 . 66 

0.458 

There  is  a  definite  decrease  in  exponent  as  NR  j^/Nr  p 
increases.  While  an  exponent  of  NR  p  that  varied  with  NR  M/NR  p 
could  have  been  used  in  the  correlation,  it  was  felt  that  this  was  an 
unnecessary  refinement.  An  average  value  of  M=  0.55  was  found 
to  correlate  the  data  within  the  experimental  uncertainty. 

The  heat  transfer  results  using  the  correlation  parameter 
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0.55 

^St^R,F  are  Presented  in  Fig.  32.  The  four  to  one  variation 
in  Nr,  p  is  seen  to  be  very  well  correlated.  Additionally,  a  least 
squares  straight  line  can  be  passed  through  all  of  the  data  (ex¬ 
cluding  runs  1-4  and  19  -  22  as  discussed  previously).  The  heat 
transfer  is  described  by  the  relation: 

NStNR,F  5  =  23.81  +  4.068  (NRjM/NR,p)  [  iU  ] 

The  standard  deviation  of  the  data  about  the  correlation  is  -  2  per 
cent.  The  95  per  cent  confidence  limits  are  t  1.3  per  cent. 

These  values  are  well  within  the  estimated  experimental  uncertainty 
of  t  4.0  per  cent  on  NstNR,F°'^5  (Appendix  A) . 

The  predicted  heat  transfer  results  (Table  III,  p.  96),  evalu¬ 
ated  at  the  extremes  of  the  NR  p  range  tested,  are  also  shown  on 
Fig.  32.  The  data  is  50  to  100  per  cent  higher  than  the  predicted. 

The  effect  of  disk  speed  (Nr  jyj-)  is  much  less  than  predicted,  as 
over  the  range  nr,m/^R(F  =  0  to  0.7,  the  predicted  heat  transfer 
increases  by  2  30  per  cent,  while  the  actual  rise  is  only  13  per  cent. 

The  actual  heat  transfer  is  surprisingly  insensitive  to  disk  speed. 

The  effect  of  adding  the  labyrinth  seal  at  the  sealing  section 
and  the  vaned  diffuser  is  significant  only  at  low  NR  m/Nr  p. 

The  earlier  data  (Runs  1  -  11  and  IS  —  2  2)  are  15  to  20  per  cent  low 
at  nR,m/Nr(f  =  0.1,  but  are  not  measurably  different  at 
^R,m/Nr,F  =  0-5.  This  reinforces  the  acceptance  of  the  present 
data  as  representative  of  rotating  disk  heat  transfer. 
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At  the  higher  Nr  j^/Nr  p,  the  fact  that  flow  does  bypass  the 
rotating  screens  has  probably  a  small  effect  on  the  heat  transfer. 
The  effect  of  the  bypass  flow  is  to  reduce  the  fluid  angular  momen¬ 
tum,  producing  lower  relative  tangential  velocities  between  the  disk 
and  the  air  stream.  This  is  analogous  to  testing  at  a  disk  speed 
(ND  A/r)  lower  than  the  actual  value,  and  the  heat  transfer  is  rela- 
tively  insensitive  to  Nr  changes  (Fig.  32).  For  the  present 
tests,  the  actual  angular  momentum  is  about  70  per  cent  of  that 
for  no  bypass  flow.  A  corresponding  30  per  cent  reduction  in  the 
test  value  of  NR  m/Nr  p  would  decrease  the  heat  transfer  by 
about  6  per  cent.  Complete  elimination  of  the  bypass  flow  would 
then  be  expected  to  increase  the  heat  transfer  coefficients  on  the 
order  of  5  to  10  per  cent  over  the  test  range. 

The  heat  transfer  coefficients  of  Smith  (19)  for  radial  inflow 
are  on  the  order  of  one-tenth  as  large  as  the  present  results.  The 
exponent  on  ND  „  is  -  0.35,  compared  to  the  present  correlation 
exponent  of  +0.55.  As  Smith  tends  to  disclaim  his  radial  inflow 
results,  a  detailed  comparison  of  the  heat  transfer  results  was  not 
made . 

In  Fig.  33,  local  heat  transfer  results  for  two  values  of 
Nr  m/Nr  p  are  plotted  as  a  function  of  radius  (r*) .  The  pre¬ 
dicted  behavior  for  NR  M/NR  p  =  0.1  is  shown  for  comparison. 

The  data  are  not  self-consistent,  and  show  little  agreement  with 
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the  magnitudes  and  trend  of  the  analysis.  These  results  are  typical 
of  the  attempts  to  obtain  local  heat  transfer  results. 


These  discrepancies  in  local  heat  transfer  behavior  are  the 
result  of  inaccuracies  in  evaluating  the  terms  in  the  series  solution 
for  the  test  surface  heat  flux.  The  term  in  Eq.  [  12  ] 
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contains  the  test  data.  If  the  second  term  in  the  integrand  is 
small  relative  to  q^”  (n)  for  the  first  few  eigenvalues  (  am  )( 
the  value  of  the  integral  will  be  small  and  the  series  will  converge 
rapidly.  This  condition  was  not  met  for  the  heat  transfer  tests,  and 
the  higher  order  terms  predominated.  The  uncertainty  in  the  numeri¬ 
cal  integration  for  these  higher  order  terms  was  large,  and  the  series 
could  not  be  accurately  summed.  The  convergence  would  be  satis¬ 
factory  and  the  local  heat  transfer  could  be  evaluated  if: 

a)  The  test  surface  temperature  variation  tQ  (  n  )  could  be 
reduced  to  +  0.2  F  from  the  present  test  variation  of 

_J_  —I 

-IF.  However,  this  reduction  was  not  feasible  with 
the  present  method  of  incident  radiant  energy  control  and 
the  present  surface  temperature  measurement  uncertainty 
of  ±  0 . 3  F . 

b)  The  heat  flux  could  be  increased  to  5000  to  7000  BTU/hr 
ft2  from  the  test  values  of  1000  to  1400  BTU/hr  ft2.  The 
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present  temperature  limitations  of  the  disk  and  the  air 
flows  obtainable  precluded  such  an  increase. 

c)  The  test  disk  thickness  could  be  reduced  to  0.1  to  0.2" 
from  the  present  thickness  of  1 . 682"  . 

The  following  conclusions  can  be  drawn  from  the  radial  inflow 
heat  transfer  results: 

1.  The  experimental  technique  for  obtaining  average  heat 
transfer  coefficients  is  satisfactory. 

2.  The  heat  transfer  did  not  correlate  as  predicted.  However, 
the  data  was  correlated  over  the  entire  test  range  by  the 
relation: 

NStNR,F°'55  =  23'81  +  4‘068  (NR,m/NR,F>  ^ 

The  spread  of  the  data  using  this  correlation  is  less 
than  the  estimated  experimental  uncertainty  of  t  4.0 
per  cent . 

3.  The  experimentally  determined  heat  transfer  is  higher 
and  less  sensitive  to  NR  m^R  F  t'nan  t^iat  predicted. 

The  possible  reasons  for  these  differences  are  discussed 
in  Section  VII. 

4.  The  effects  of  the  test  section  inlet  and  bypass  flow  are 
small,  and  confined  to  low  NR  M/NR  p . 
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5.  Local  heat  transfer  coefficients  could  not  be  satisfac¬ 
torily  obtained.  An  increase  in  test  heat  flux  or  a 
decrease  in  disk  thickness  would  be  necessary  to  obtain 
local  heat  transfer  behavior. 


FIGURE  30 


AVERAGE  KEAT  TRANSFER  RESULTS  FOR 
RADIAL  INFLOW 

0.2 

NqtND  „  AS  A  FUNCTION  OF  N  /N  „ 
ft ,  r  R.MR,r 


FIGURE  31 

AVERAGE  HEAT  TRANSFER  RESULTS  FOR 
RADIAL  INFLOW 

N<;t  AS  A  FUNCTION  OF  N 
bt  R,F 

(a)  Nr  m/Nr  f  =  0.095,  0.31,  0.66,  and  1.27 

(b)  NR  p^/N r  P  =  0.15,  0.50,  0.99,  and  1.66 
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FIGURE  32 


AVERAGE  HEAT  TRANSFER  RESULTS  FOR 


RADIAL  INFLOW 


NstNR  f°'55  AS  A  function  of  nr  m/nr,f 
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FIG. 32 


FIGURE  33 


LOCAL  HEAT  TRANSFER  RESULTS 
FOR  RADIAL  INFLOW 

0 . 2 

Nst  inr.f  as  a  FUNCTION  OF  r* 
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TABLE  I 

RADIAL  INFLOW  HEAT  TRANSFER  OATA 
*o  -  0. 90" 


RUN 

W  p 

Ib/hr  lb/ft3 

(0 

RPM 

C«ir 

F 

*dUk 

F 

h 

BTU 

hr  ft2 

N  _ 
R.F 

“  xlO" 

nr,m 

'5  xlO' 

nr,f 

5 

19 

1691 

0.0736 

174 

76.  1 

125.3 

20.16 

5.40 

0. 479 

0.069 

20 

1660 

0.0730 

140 

60.8 

134.8 

20.92 

5.33 

0.380 

0.071 

21 

1982 

0.0729 

220 

76.  1 

133.9 

21.91 

6.31 

0.600 

0.095 

22 

1974 

0.0728 

216 

79.4 

135.5 

22.04 

6.27 

0.592 

0.094 

1 

2099 

0.0721 

275 

63.1 

139.7 

20.77 

6.64 

0.736 

0.111 

2 

2091 

0.0720 

280 

64.3 

140.1 

20.70 

6.60 

0.747 

0.113 

207 

796 

0.0736 

140 

77.4 

126.4 

16.74 

2.36 

0.365 

0. 164 

208 

795 

0.0734 

143 

78.8 

129.0 

17.00 

2.35 

0.392 

0.  167 

216 

1469 

0.0732 

220 

77.6 

118.3 

21.43 

4.41 

0.602 

0. 137 

219 

1465 

0.0730 

225 

78.7 

121.7 

21.47 

4.39 

0.613 

0.140 

220 

1872 

0.0725 

324 

80.6 

125.7 

23.98 

5.51 

0.674 

0. 159 

221 

1666 

0.0723 

330 

62.3 

127.4 

23.66 

5.49 

0.886 

0. 162 

3 

2110 

0.0730 

332 

78.4 

123.6 

22.00 

6.72 

0.905 

0.135 

4 

2090 

0.0726 

345 

60.6 

126.5 

21.59 

6,63 

0.932 

0.141 

170 

623 

0.0726 

196 

83.4 

136.5 

15.47 

1.629 

0.526 

0.289 

169 

623 

0.0726 

194 

83.4 

138.4 

15.50 

1.826 

0.523 

0.266 

166 

915 

0.0726 

352 

83.1 

136.4 

16.26 

2.69 

0.946 

0.353 

167 

916 

0.0726 

349 

62.8 

137.8 

16.39 

2.69 

0.941 

0.350 

166 

1250 

0.0727 

396 

60.6 

136.2 

20.83 

3.66 

1.072 

0.291 

165 

1251 

0.0729 

394 

79.4 

134.9 

20.79 

3.69 

1.071 

0.290 

1« 

1655 

0.0730 

549 

77.4 

129.6 

22.77 

4.90 

1.498 

0.306 

163 

1656 

C.0731 

550 

76.3 

126.1 

23.00 

4.91 

1.506 

0.307 

161 

1871 

0.0726 

645 

78.0 

129.3 

24.50 

5.53 

1.750 

0.316 

160 

1873 

0.0727 

640 

77.0 

129.0 

24.13 

5.55 

1.742 

0.314 

209 

773 

0.0730 

458 

82.0 

133.3 

17.38 

2.27 

1.243 

0.547 

210 

827 

0.0726 

449 

63.0 

133.3 

17.72 

2.43 

1.214 

0.500 

211 

1130 

0.0725 

642 

84.2 

131.5 

20.43 

3.31 

1.726 

0.521 

212 

1126 

0.0725 

654 

64.9 

132.8 

20.18 

3.31 

1.754 

0.531 

222 

1660 

0.0720 

1025 

84.4 

130.5 

24.68 

5.45 

2.73 

0.501 

223 

1858 

0.0719 

1026 

85.1 

130.2 

25.21 

5.44 

2.73 

0.502 

9 

1795 

0.0729 

1008 

78.0 

124.0 

25.52 

5.71 

2.75 

0.460 

10 

1789 

0.0729 

1010 

78.0 

125.8 

24.72 

5.70 

2.75 

0.483 

11 

1790 

0.0724 

1005 

81.3 

131.8 

25.12 

5.67 

2.71 

0.477 

172 

589 

0.0735 

389 

78.0 

137.7 

15.25 

1.741 

1.069 

0.614 

173 

587 

0.0734 

392 

79.0 

138.5 

15.35 

1.734 

1.074 

0.619 

174 

729 

0.0732 

505 

80.2 

137.6 

16.65 

2.15 

1.377 

0.641 

175 

728 

0.0731 

507 

80.9 

137.7 

17.04 

2.14 

1.379 

0.643 

176 

1123 

0.0736 

852 

77.0 

130.6 

21.12 

3.33 

2.34 

0.705 

177 

1122 

0.0734 

872 

77.9 

131.2 

21.16 

3.32 

2.39 

0.721 

178 

1519 

0.0727 

1078 

81.6 

134.9 

23.79 

4.47 

2.92 

0.652 

179 

1516 

0.0726 

1083 

62.6 

136.4 

23.68 

4.45 

2.92 

0.655 

5 

1611 

0.0736 

1310 

77.0 

131.4 

25.99 

5.15 

3.62 

0.703 

6 

1594 

0.0731 

1295 

78.  1 

132.9 

25.78 

5.07 

3.53 

0.697 

7 

1604 

0.0729 

1275 

79.8 

134.4 

25.60 

5.09 

3.46 

0.660 

160 

1862 

0.0720 

1282 

85.0 

133.6 

25.63 

5.45 

3.42 

0.627 

161 

1859 

0.0719 

1275 

65.8 

134.3 

25.80 

5.44 

3.39 

0.623 

171 

609 

0.0737 

528 

63.4 

138.5 

16.00 

1.605 

1.457 

0.807 

213 

563 

0.0727 

598 

85.2 

139.4 

16.42 

1.648 

1.606 

0.975 

214 

563 

0.0727 

599 

85.2 

140.4 

16.12 

1.650 

1.610 

0.976 

215 

994 

0.0727 

1087 

83.8 

132.9 

20.86 

2.92 

2.93 

1.005 

224 

1744 

0.0719 

1912 

65.7 

127.4 

27.67 

5.10 

5.06 

0.996 

225 

1744 

0.0719 

1915 

85.2 

128.6 

26.63 

5.11 

5.10 

0.998 

190 

514 

0.0733 

705 

79.9 

134.6 

15.77 

1.515 

1.926 

1.271 

189 

515 

0.0733 

668 

79.6 

135.0 

15.63 

1.519 

1.660 

1.238 

166 

863 

0.0737 

1170 

75.8 

127.5 

20.32 

2.56 

3.23 

1.262 

167 

666 

0.0739 

1163 

74.6 

126.5 

20.40 

2.57 

3.23 

1.255 

166 

1043 

0.0728 

1518 

80.9 

134.7 

22.63 

3.07 

4.12 

1.339 

163 

1363 

0.0730 

2000 

78.3 

127.9 

25.43 

4.05 

5.45 

1.346 

162 

1736 

0.0733 

2242 

74.3 

120.7 

27.59 

5.22 

6.17 

1.163 

191 

634 

0.0731 

1125 

81.1 

136.5 

17.81 

1.666 

3.06 

1.636 

192 

634 

0.0731 

1132 

80.9 

136.3 

17.61 

1.867 

3.08 

1.649 

216 

945 

0.0728 

1710 

82.3 

130.6 

22.42 

2.61 

4.62 

1.646 

217 

955 

0.0728 

1696 

62.2 

130.9 

22.23 

2.61 

4.59 

1.636 

6 

1064 

0.0727 

2147 

61.4 

131.8 

24.46 

3.37 

5.61 

1.721 

0.01462 
0.01526 
0.01356 
0.01369 
0.01212 
0.012  14 

0.0277 

0.0262 

0.01899 

0.01905 

0.01669 

0.01666 

0.01276 

0.01270 

0.0328 

0.0328 

0.0263 

0.0265 

0.0220 

0.0219 

0.01615 

0.01630 

0.01726 

0.01700 

0.0297 

0.0263 

0.0236 

0.0236 

0.01749 

0.01789 

0.01744 

0.01694 

0.01721 

0.0342 

0.0345 

0.0305 

0.0308 

0.0246 

0.0249 

0.0207 

0.0206 

0.01941 

0.01980 

0.01957 

0.01830 

0.01830 

0.0347 

0.0385 

0.0378 

0.0277 

0.0211 

0.0201 

0.0405 

0.0400 

0.0311 

0.0311 

0.0266 

0.0245 

0.0207 

0.0370 

0.0370 

0.0310 

0.0307 

0.0262 


0.205 
0.214 
0. 1960 
0. 1975 
0.1771 
0. 1770 

0.329 
0.335 
0.255 
0.256 
0.237 
0.237 
0. 1872 
0. 1869 

0.369 

0.369 

0.320 

0.322 

0.266 

0.264 

0.250 

0.252 

0.243 

0.239 

0.352 

0.338 

0.302 

0.300 

0.246 

0.252 

0.247 

0.240 

0.243 

0.382 

0.385 

0.355 

0.360 

0.316 

0.318 

0.279 

0.276 

0.270 

0.274 

0.271 

0.257 

0.257 

0.324 

0.425 

0.418 

0.343 

0.292 

0.279 

0.440 

0.435 

0.375 

0.375 

0.356 

0.324 

0.268 

0. 420 
0.420 
0.361 
0.378 
0.360 


20.8 

21.5 
21.0 
21.1 
19.31 
19.30 

25.0 

25.4 

24.2 

24.2 

24.4 

24.2 

20.6 

20.2 

25.7 

25.7 

25.5 

25.7 

25.4 

25.3 

24.5 

24.7 

24.4 
24.0 

26.2 

26.0 

25.8 

25.6 

24.8 

25.4 

25.6 

24.8 

25.2 

26.1 

26.3 

26.1 

26.4 
27.0 

27.1 

26.5 

26.4 

26.8 

27.2 
27.0 

26.2 

26.2 

27.0 

28.5 
26.0 
26. 1 
29.1 

27.8 

26.6 

26.3 

29.4 

29.4 

29.8 

29.6 

26.8 

29.4 

29.4 

30.8 

30.5 

30.9 
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C .  Radial  Outflow  Heat  Transfer  Results 

At  the  conclusion  of  the  radial  inflow  studies,  a  brief  investi¬ 
gation  of  the  heat  transfer  for  forced  radial  outflow  was  conducted. 
The  purpose  of  this  study  was  to  obtain  heat  transfer  coefficients 
and  compare  them  with  the  recently  available  data  of  Kreith, 
Doughman,  and  Kozlowski  (10).  A  correlation  obtained  by  Smith 
(IS)  also  affords  a  comparison.  These  comparisons  help  in  the  test 
facility  evaluation. 

To  obtain  a  radially  outward  flow  of  air  through  the  test 
section,  the  blower  connections  were  reversed.  The  blower  induced 
air  from  the  atmosphere  and  discharged  into  the  metering  section. 

The  vaned  diffuser  was  removed  as  it  would  have  imparted  a  swirl 
to  the  air  entering  the  test  section.  No  other  modifications  were 
made . 

Heat  transfer  coefficients  were  obtained  at  a  single  shroud- 
to-disk  clearance  ratio  of  zQ  =  0.113  over  a  flow  Reynolds  number 
range  NR  p  =  1.5  x  105  to  6.5  x  105  for  Reynolds  number  ratios 
of  Nr  M/NR  F  =  °-30'  0-58,  and  1.14.  The  average  heat  transfer 
coefficient  was  determined  in  the  same  manner  as  for  radial  inflow. 
To  afford  a  direct  comparison  with  the  results  of  (10),  the  heat 
transfer  was  characterized  by  the  average  Nusselt  number.  The 


h _ Tq 

k 


q  /  ^  1  disk 
k 
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The  average  heat  transfer  data  and  results  are  tabulated  in 
Table  II.  The  results  are  grouped  according  to  the  nominal  test 
value  of  Nr  m/Nr  f‘  The  correlation  parameter  suggested  in 
(10)  is  designated  as  K. 

In  Fig.  34,  the  average  Nusselt  number,  NNu,  is  presented 

as  a  function  ot  ND  „.  The  equation  of  the  least  squares  straight 
K,  Jr 

line  passing  through  all  of  the  data  is 

NNu  =  0,312  nR,F 

The  confidence  limit  of  the  data  about  this  line  is  ±  3.5  per  cent. 
This  is  about  equal  to  the  estimated  experimental  uncertainty  in  NNu  , 
of  -  3.4  per  cent.  While  there  does  appear  to  be  a  slight  decrease 
in  with  increased  NR  M/NR  p<  the  magnitude  of  the  decrease 

is  about  equal  to  the  experimental  uncertainty. 

In  Fig.  13  of  (10),  Sherwood  number  (the  mass  transfer  para¬ 
meter  corresponding  to  the  Nusselt  number)  is  plotted  as  a  function 

of  Nn  )  with  N  as  a  parameter.  The  shroud-to-disk 

h,  r  R„ M 

spacing  is  zQ*  =  0.025.  These  results  were  crossplotted  to 
obtain  Nj^u  as  a  function  of  NR  p  for  the  present  nominal  test 
values  of  Np  M/NR  p,  and  are  also  shown  in  Fig.  34. 

The  present  results  are  seen  to  be  10  to  50  per  cent  lower  than 
those  of  (10)  over  the  test  range.  Tne  slope  of  the  present  correla¬ 
tion  is  0.626,  while  those  from  (10)  are  0.76,  0.84,  and  1.05  for 
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Nr  m/Nr  p  =  0.30,  0.58,  and  1.14,  respectively.  The  present 
heat  transfer  results  are  less  sensitive  to  Nr  p  changes  than  those 
of  (10) . 

It  is  surprising  that  there  is  no  effect  of  disk  rotation  (Np  jyj) 
on  the  present  heat  transfer  results.  Increased  disk  speeds  would 
be  expected  to  increase  the  heat  transfer.  The  results  of  (10)  show 
that  the  heat  transfer  increases  25  to  50  per  cent  over  the  Nr  ^ 
range  tested . 

Over  the  zq*  range  (0.012  to  0.02  5)  of  (10),  the  heat  trans- 

*  a  r  r 

fer  is  correlated  by  zQ  u  .  Extending  this  correlation  to  the 
present  data  where  zQ*  =  0.113  indicated  that  the  present  heat 
transfer  is  still  2  5  to  80  per  cent  low.  These  discrepancies  must 
be  resolved  by  further  tests. 

The  radial  outflow  correlation  obtained  by  Smith  (19)  was 

evaluated  for  the  present  test  conditions  and  is  also  presented  in 

Fig.  34.  The  present  data  is  again  seen  to  be  low.  The  agree- 

,  .  0.75 

ment  in  slope  is  better;  the  heat  transfer  of  (19)  varies  as  Np^p 
while  the  present  variation  is  N  0.626.  A  detailed  comparison 
is  not  too  meaningful,  however,  due  to  the  large  uncertainties 
surrounding  the  test  procedure  of  Smith  (19). 

In  (10),  all  of  the  results  are  correlated  by  the  relation: 
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N 


Sherwood 


(  =  NNu) 


0.55 


1.36  +1.29 


'nR,M 

105 


+  3.571^$  .  3.51  1^0  + 
(WR,F  ,(0-83-  0.12  ^M_) 


[IT] 


The  results  are  presented  as  N  as  a  function  of  K  in  Fia  35 

Nu  ‘  ' 

and  are  again  seen  to  be  50  to  70  per  cent  low.  This  is  compatible 

with  the  previous  data  comparisons. 

The  heat  transfer  behavior  using  the  correlation  parameter  K 

of  Eq.  [  17  ]  is  rather  peculiar.  This  was  seen  by  evaluating  NNu 

using  the  present  correlation,  Eq.  [  16],  for  NR  M/NR  p  =  1. 

This  result  is  also  shown  on  Fig.  35,  and  it  is  seen  that  for  a 

given  value  of  K,  there  are  two  possible  values  of  N  .  While 

Nu 

the  present  results  are  at  some  variance  with  those  of  (10),  they 
do  show  the  same  trends  on  and  NR  R  coordinates.  These 

discrepancies  are  magnified  using  the  correlation  parameter  of 
Eq-  [  17  ]  .  This  reduces  the  usefulness  of  such  a  correlation 
for  comparison  purposes. 

The  possible  reasons  for  the  discrepancy  between  the  present 
results  and  those  of  Kreith,  et  al  (10)  are: 

a)  The  present  zq*  is  about  four  times  the  maximum 

tested  in  (10).  The  effect  of  zD*  as  determined  in  (10) 
may  not  extend  outside  the  test  range. 
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b)  The  inlet  conditions  may  differ.  In  the  present  test  flow 
facility,  the  air  flow  passage  had  a  right  angle  bend  six 
diameters  upstream  of  the  test  section  inlet.  This  could 
affect  the  flow  pattern  entering  the  test  section. 

c)  The  flow  pattern  in  the  present  test  section  is  undoubtedly 
different  from  that  of  (10)  due  to  the  presence  of  the 
rotating  screens  and  sealing  section.  The  screens  proba¬ 
bly  reduce  the  air  velocity  next  to  the  disk,  while  the  by¬ 
pass  flow  in  the  sealing  section  probably  increases  the 
air  velocity  next  to  the  shroud.  As  a  result,  the  heat 
transfer  would  be  low  as  indicated. 

The  conclusions  that  can  be  drawn  from  the  radial  outflow  heat 
transfer  investigations  are: 

1.  The  present  results  show  no  effect  of  disk  rotation  (NR  M) 

on  the  heat  transfer.  The  heat  transfer  is  correlated  by 

0.626 

NNu  *  °'312  NR,F  tl6] 

2.  The  heat  transfer  coefficients  are  25-80  per  cent  lower 

than  those  of  (10) . 

3.  To  accurately  establish  the  radial  outflow  heat  transfer 

behavior,  additional  results  must  be  obtained.  The 
effects  of  the  inlet,  screens,  and  sealing  section  must 
be  determined . 
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FIGURE  34 


COMPARISON  OF  AVERAGE  HEAT  TRANSFER 
RESULTS  FOR  RADIAL  OUTFLOW 

NNu  AS  a  FUNCTION  OF  Nr  f 


FIGURE  35 

AVERAGE  HEAT  TRANSFER  RESULTS  FOR 
RADIAL  OUTFLOW 

NNu  AS  A  FUNCTION  OF  THE  CORRE¬ 
LATION  PARAMETER,  K 
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FIG.  35 


TABLE  II 


RADIAL  OUTFLOW  HEAT  TRANSFER  DATA 


z  -  0.90” 
o 


RUN 

V 

P 

CD 

tair 

Cdl«k 

h 

Vf 

nr,m 

K 

lb/hr 

lb/ft3 

RFM 

F 

F 

BTU 

xlO'5 

xlO'5 

^.F 

hr  ft2  F 

236 

609 

0.0712 

201 

97.2 

147.2 

14.85 

1.756 

0.521 

0.297 

636 

1420 

237 

609 

0.0711 

203 

98.0 

148.3 

14.76 

1.753 

0.525 

0.299 

632 

1417 

238 

873 

0.0709 

296 

99.9 

142.5 

17.86 

2.51 

0.762 

0.304 

762 

1814 

239 

883 

0.0709 

294 

100.0 

142.8 

17.77 

2.54 

0.756 

0.298 

757 

1831 

232 

1331 

0.0725 

440 

88.9 

130.1 

23.53 

3.88 

1.175 

0.303 

1019 

2336 

233 

1314 

0.0722 

442 

91.3 

131.4 

24.19 

3.82 

1.171 

0.307 

1042 

2311 

226 

1917 

0.0712 

620 

101.3 

134.5 

29.32 

5.49 

1.600 

0.291 

1249 

3028 

227 

1913 

0.0712 

620 

101.8 

135.5 

28.90 

5.48 

1.598 

0.292 

1230 

3023 

240 

427 

0.0709 

263 

98.6 

152.4 

11.73 

1.228 

0.678 

0.552 

501 

1076 

241 

427 

0.0710 

258 

98.2 

151.8 

11.56 

1.228 

0.666 

0.542 

494 

1076 

242 

937 

0.0707 

615 

88.5 

135.5 

18.53 

2.69 

1.580 

0.586 

791 

1911 

243 

932 

0.0704 

626 

91.2 

138.2 

18.53 

2.67 

1.594 

0.597 

788 

1908 

!  234 

1448 

0.0719 

96  5 

93.8 

134.4 

23.24 

4.19 

2.54 

0.606 

1000 

3321 

235 

1444 

0.0716 

938 

96.1 

135.8 

23.78 

4.17 

2.45 

0.588 

1018 

3247 

228 

1930 

0.0711 

1240 

102.3 

134.4 

30.09 

5.52 

3.19 

0.577 

1280 

4002 

229 

1929 

0.0710 

1240 

102.6 

134.0 

30.77 

5.52 

3.18 

0.577 

1309 

4004 

244 

435 

0.0694 

577 

91.8 

151.0 

11.11 

1.233 

1.435 

1.164 

468 

1108 

245 

434 

0.0693 

584 

92.8 

151.0 

11.07 

1.229 

1.450 

1.180 

466 

1110 

246 

717 

0.0698 

913 

94.5 

144.2 

14.45 

2.04 

2.29 

1.124 

611 

2094 

247 

716 

0.0698 

920 

94.7 

143.5 

14.71 

2.04 

2.31 

1.133 

622 

2109 

249 

1086 

0.0699 

1467 

97.9 

139.9 

18.35 

3.09 

3.69 

1.193 

776 

2993 

230 

1727 

0.0705 

2152 

105.0 

138.9 

25.03 

4.93 

5.47 

1.109 

1060 

1748 

231 

1739 

0.0710 

2147 

100.9 

134.1 

26.12 

4.99 

5.52 

1.107 

1)12 

1702 

K  Is  the  correlation  parameter  of  (10). 
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VI  ANALYTICAL  INVESTIGATION 


The  specific  objective  of  the  analytical  investigation  is  to  pre¬ 
dict  the  heat  transfer  behavior.  As  achievement  of  such  a  goal  is 
optimistic,  it  is  hoped  that  at  least  the  study  will  generate  the  param¬ 
eters  governing  the  he;at  transfer,  indicate  the  expected  trends,  and 
allow  some  insight  into  the  experimentally  determined  fluid  dynamics 
and  heat  transfer  behavior. 

Integral  methods  will  be  used  to  solve  the  differential  equations 
governing  the  fluid  dynamics  and  heat  transfer.  The  Reynolds  analogy 
will  be  used  to  relate  the  transfer  of  heat  to  the  transfer  of  momentum. 
Boundary  layers  will  be  assumed  to  exist  on  the  disk  and  shroud.  The 
procedure  will  be  to  first  solve  for  the  fluid  behavior  external  to  the 
boundary  layers .  The  momentum  integral  equation  will  then  be  solved 
to  yield  the  hydrodynamic  boundary  layer  growth  and  the  wall  shear 
stress.  The  Reynolds  analogy  and  the  energy  integral  equation  will 
then  be  solved  to  evaluate  the  heat  transfer  from  the  rotating  disk. 

The  following  assumptions  are  made  in  the  analysis: 

1.  The  flow  and  heat  transfer  patterns  are  axisymetric  and 

steady;  there  are  no  property  variations  in  the  tangential 
direction . 

2.  The  fluid  behaves  as  an  uncompressed  gas;  density  changes 

due  to  pressure  and  temperature  changes  are  small. 

3  .  The  flow  is  turbulent . 


78 


4.  The  boundary  layer  thickness  is  small  compared  to  the  disk- 

to- shroud  spacing.  The  boundary  layers  are  thin  enough  so 
that  the  fluid  pressure  varies  with  radius  only. 

5.  The  fluid  velocity  external  to  the  boundary  layers  is  not 

influenced  by  the  boundary  layers  (essentially  zQ*  =  00  ). 

6.  The  wall  shear  stress  is  given  by  the  Blasius  relation.  This 

relation  has  been  found  to  adequately  describe  the  wall 
shear  stress  for  many  boundary  layer  flows. 

7.  The  velocity  and  temperature  profiles  are  described  by  l/7th 

power  profiles.  The  l/7th  power  profile  has  been  found  to 
adequately  describe  the  boundary  layer  profile  in  many 
boundary  layer  flows,  and  is  compatible  with  the  Blasius 
relation  for  wall  shear  stress. 

8.  The  eddy  diffusivities  for  heat  {  eH  )  and  momentum  (em) 

are  equal  over  the  boundary  layer.  This  implies  that 
Npr=  1. 

S.  The  fluid  temperature  at  the  edge  of  the  thermal  boundary 
layer  and  the  disk  surface  temperature  do  not  vary  with 
radius . 

These  assumptions  are  similar  to  those  usually  made  in  studies 
of  boundary  layer  flows.  The  most  questionable  ones  are  those  con¬ 
cerning  the  shear  stress  relations  and  the  velocity  and  temperature 
profiles.  The  present  assumptions  are  made  as  no  information  is 
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available  on  the  profiles  for  radial  inflow  over  a  rotating  disk. 

In  an  integral  analysis,  the  heat  transfer  is  not  too  sensitive  to 
the  choice  of  profile  as  the  profiles  are  integrated  over  the  boundary 
layers.  The  solution  is  more  sensitive  to  the  shear  stress  relation 
and  external  flow  pattern.  The  solution  and  assumptions  will  have 
to  be  verified  experimentally. 


A.  Behavior  of  the  Fluid  External  to  the  Boundary  Layers 

The  fluid  behavior  external  to  the  boundary  layers  on  the  disk 
and  shroud  is  described  by  the  following  equations  derived  for  a 
stationary  control  volume  located  in  the  fluid.  The  coordinate 
system  is  described  in  Fig.  1,  p.  10. 

Geometry: 


d  8 


v 

u 


dr 


Continuity  equation: 

1  3  (ru) 

r  3  r 


-  0 


[18a] 


[19] 


Momentum  equation  in  the  radial  direction: 

-L_(  u21L  -  Z?)  =  _  _L 

gc  ar  r  P  3'r 

Momentum  equation  in  the  tangential  direction: 
1  u  3  (rv)  _ 


[20] 


[21] 
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Boundary  conditions: 


At  r  =  ro  u(rQ,  z) 

v(ro<  2) 


W 


2  r  z  p 
o  o 


u 


At  z  =  0  u  (r,  0)  =0 

v  (r,  0)  =  r  u 

The  boundary  conditions  correspond  to  those  both  for  radial 
inflow  turbines  and  for  the  present  model.  The  relative  fluid  velocity 
is  radial  with  respect  to  the  disk  at  the  disk  periphery. 

Eqs .  [19]  and  [21]  can  be  integrated  directly.  Using  the  boundary 

conditions,  the  fluid  velocities  are  written  as: 


1  ,  -  w 

r  2  n  zQ  p 

-)  -  ->£(«,  1 

2  IT  ro  r* 

[22] 

1  ,  2  . 

~  —  (  "o  “  }  = 

NR,M  (rQ  }  r* 

[23] 

The  tangential  velocity  component,  v,  exhibits  what  is  called  free 
vortex  flow.  The  angle  of  the  resultant  velocity  with  respect  to 
the  radial  direction  is 


tan  (  a  )  =  — 

u 


N 


2  TI 


R,M 


N 


R,F 


[2U] 


These  results,  [22]  ,  [23],  and  [  24  ] ,  are  presented  in  Figs . 

14  to  28  in  comparison  with  the  experimental  results.  The  pressure 
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gradient  is  obtained  from  L  20  using  the  fluid  velocities  [22] 
and  [23]: 


3c 


dp 

dr 


Np  -p  \  2  . ; 

-r'r  +  <nr,mI 


or: 


gc  dP 

o  r  dr* 
o 


(r!r")2  (1+p2 


[25] 


Eq.  [  25  j  is  presented  in  Fig.  2S  in  comparison  with  the  experi¬ 
mentally  determined  wall  static  pressure  variation.  Eq.[  18  ]  can 
now  be  solved  to  compute  the  path  of  the  fluid  as  it  flows  radially 
inward.  It  is  more  illuminating,  though, to  study  the  fluid  path 
with  respect  to  the  disk.  The  relative  tangential  velocity  and 
relative  flow  angle  are  related  by: 

da  =  [18b] 

R  u  r 


where 


r  ... 


=  N 


(1  -  r*  2  ) 


[26] 


R,M 


The  relative  angle  is  then  given  by: 


2 1.  NR,M 


N 


-  1 


+  In 


[27] 


R,F 


L 


Eq.  [  27  ]  is  plotted  in  Fig.  36  for  various  values  of  p. 

For  Nr  m/Nr  f  <  0.1,  there  is  little  effect  of  disk  rotation  on 


the  fluid  path . 
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B.  Hydrodynamic  Solution 


The  velocities  and  pressure  gradient  obtained  for  the  fluid 

external  to  the  boundary  layer  can  now  be  used  to  help  solve  the 

hydrodynamic  problem  inside  of  the  boundary  layer .  The  governing 

equation  is  the  radial  momentum  differential  equation  integrated 

over  the  hydrodynamic  boundary  layer  (the  integration  is  carried  out 

in  Appendix  F) .  This  equation  is  written  for  a  stationary  control 

volume  located  in  the  fluid:  * 

6 


where  the  coordinates  and  the  boundary  layer  thickness  are 
normalized  by  the  disk  radius. 

The  pressure  gradient  is  given  by  [  25]  .  The  fluid  veloci¬ 
ties  inside  of  the  boundary  layer  are  assumed  to  vary  from  their 
absolute  value  at  the  disk  to  their  absolute  value  at  the  edge  of 
the  boundary  layer  with  the  l/7th  power  of  distance: 


u  = 


i¥  <— >  - 

2  *  rQ  r 


1/7 

) 


v  = 


N 


R,M 


(1  -  r*2) 
r* 


[29] 


[30] 


and  the  relative  tangential  velocity  is 


83 


84 


85 


The  differential  equation  for  6  *  can  be  obtained  by  substi¬ 


tuting  [  25  ],[  34]  ,  [35]  [36],  and  [37]  into  [28].  After 


multiplying  by  r*  and  dividing  by 


becomes: 


7  i/4  1/4 

0.0225  (4rr-=)  '  r 


NR,F 


N 


R,F 


2  n 


1/4  L_ 


1  +  e 2  (1  -  r*2  )2 


,  Eq.  [  28: 


n  3/8 


(  1  +  6  )  = 


d6 


8 


dr* 


-  8 


dr* 


+  7r*A  +  28) 
36 


[38a] 


Combining  terms, 

7  d  6  * 
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dr* 


R,F 


2  w  1/4  r*l/4 

=  0.022  5  (^— _) 


6  * 


1/4 


2  *2 


(1+s2>  ■  ~hr  8~r* 


-  _8  J 
36 


1  + 


(1  -  r*2) 


> 


3/8 


[38b] 


1/4 

Multiplying  [  38b]  by  6  *  ,  and  recognizing  that 


4/4 


_d6_ 


dr* 


d  <s  * 
dr* 


5/4 


allows  [38b]  to  be  written  as 
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dr* 


(  5*  NR,F1/5)5/4  +  (fi*  NR,F1/5)5/4 


20  (1  +82)  +  5  g  2  r*  +  _5_ 


2 

1/4 


14 


2  *3 
6  r* 


0.4579  r*  ^1+6(1 


*2  x2 

r*  ) 


2~1  3/8 


Eq.  [39]  is  a  first  order  differential  equation  linear  in  (6  *  N  1//5)5/4^ 

R  t  F 

It  is  solved  with  the  aid  of  the  integrating  factor: 


20  (1  +8  2) 
7  r* 


+  -^f-  e  r 


*  +  5 


2  _*3 


14 


dr* 


20 


U+A  > 


=  r* 


2.  2 
S  (r*  + 


14 


Eq.  [39]  can  then  be  directly  integrated. 

(  6  *  Nr(F  ^  )5/"4  is  that  at  r  =  rQ  (r*  =  1) , 


The  initial  condition  on 
1/5 

6  =  0  (6  *  N  =  0) . 

R„  F 


The  hydrodynamic  solution  is: 


fc*NR  f1/5)5/4 


0.4579 


r 


-*  ,13.  +  20  P  * 

'28  7  > 


c  2  2  4 

I  (r4+-;i-  \ 

4  v  14  ' 


(1  +  e  2)  _5_a  2  (r.2  +  -Ed  , 

e  4 


L.+  .2a-,V] 


3/8 


ihis  is  an  explicit  relation  for  6*  as  a  function  of  r*,  with  ND  t? 

K,  1 

and  B  =  2  w  NR M/NR, F  as  parameters  .  In  general,  [40]  can 
not  be  integrated  analytically,  and  a  numerical  solution  must  be 
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obtained.  For  the  case  of  no  rotation  (  B  =  2tt  p  =  0), 

[  40  ]  can  be  directly  integrated: 


1/5 

6  *  NR,F 


0.367  r* 


45 

28 


[1^1] 


In  Fig.  37,  the  hydrodynamic  boundary  layer  is  shown  as  a 
function  of  radius  for  various  values  of  p-  The  results 

for  Nr  m;/^r  f  greater  than  zero  were  obtained  using  the  Burroughs 
220  digital  computer  at  the  Stanford  Computation  Center.**  The 
boundary  layer  thickens  conventionally  at  first,  then  reaches  a 
maximum  and  decreases.  This  unusual  decrease  in  thickness  is 
due  to  fluid  acceleration.  As  the  fluid  flows  radially  inward,  the 
resultant  velocity  increases.  The  radial  component  increases  due 
to  area  reduction  and  the  tangential  c  omponent  to  the  free  vortex 
effect.  This  increasing  velocity  thins  the  boundary  layer.  The 
boundary  layer  is  thinner  at  the  higher  values  of  F  as 

the  tangential  velocity  component  is  higher. 

It  was  not  possible  to  obtain  results  for  NR  F  9reater 

than  0 . 7  on  the  computer.  This  resulted  from  the  values  of 

_  _20_  (i  +  B2) 

r*  exceeding  the  computer  capacity  at  high  B 

and  low  r*  . 
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C.  Heat  Transfer  Solution 


The  heat  transfer  from  the  rotating  disk  is  obtained  using  the 
energy  integral  equation  and  the  Reynolds  analogy  relations  .  The 
thermal  boundary  layer  thickness  is  evaluated  first.  The  energy 
differential  equation  integrated  over  the  thermal  boundary  layer  is 
(Appendix  F): 


q 


w 


P  cp 


[1*2] 


In  this  equation,  qw"  and  a  *  are  unknown.  Another  relation  is 
obtained  by  assuming  mechanistic  descriptions  for  the  transfer  of 
heat  and  momentum: 


_2L_ 


p  Cp 


a  t 

3  z 


gc  T 
p 


eM 


3 

7 


VR 

z 


'M 


.  /  2  ,  T 

l 7  U  +  VR 

3  z 


[1*3] 

[hk] 


e  H  anc*  £M  are  diffusivities  for  heat  and  momentum  in 
turbulent  flow.  The  shear  stress  relation  [  44]  ,  like  the  Blasius 
relation,  describes  the  shear  in  terms  of  the  velocities  with  respect 
to  the  disk  (relative  velocities) . 

Because  of  the  axial  symmetry  of  the  problem,  there  is  no  heat 
convected  from  the  disk  due  to  flow  in  the  tangential  direction. 
Hence  the  analogy  between  heat  and  momentum  transfer  must  be 
written  in  terms  of  the  radial  shear  stress.  The  radial  component 
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of  the  shear  stress  is 


£c  T  r  = 
P 


g 


r 


P 

(both  u  and 


3  Vr 

-  X  cos 

M  3  z 


3  u 

E  M  3  z 

vary  as 


angle  between  u 
and 


1/7 

) 


). 


[U5] 


The  Reynolds  analogy  assumes  that  the  diffusivity  for  heat 
(e^  )  equals  the  diffusivity  for  momentum  )  over  the  boundary 

layer.  The  diffusivities  can  then  be  eliminated  from  [43]  and  [44], 
and  the  wall  heat  flux  can  be  written  in  terms  of  the  wall  shear 
stress: 


,  Sc Trw. 
P 


3  t 
3  z 


3  u 
3  Z 


z  =  0 


z  =  0 


[1*6] 


Note  that  [46]  specifies  that  near  the  wall  the  temperature 
profile  must  vary  with  the  same  power  of  distance  as  the  velocity 
profile.  The  temperature  profile  will  be  assumed  to  be: 


hs.  ll.  )  = 


1/7 


[1*7] 


where  the  disk  temperature,  tw,  and  the  external  fluid  temperature, 
t  £  ,  do  not  vary  with  radius. 

There  are  now  two  relations,  [  42]  and  [  46],  for  the  wall  heat 
flux  and  the  thermal  boundary  layer  thickness.  Eq.  [46]  is 
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evaluated  first . 


at  the  wall  is: 


3  t 
3  2 
3  u 
3  z 


The  ratio  of  the  temperature  and  velocity  gradients 


[1+8] 


With  [48]  and  the  Blasius  relation  for  shear  [34],  [46]  becomes: 

2  -]  3/8 


PCP 


=  0 


r*3/4  6*i/4  q  1/7 


[1+93 


Using  the  assumed  velocity  and  temperature  profiles,  [42]  is 
integrated: 


w 


p  Cp 


=  .  (JwjJa }  L^iL)  (4U  7 

r*  i  2  »  /  fo  72  dr* 


d  (6*  -8A 


[50] 


Eqs .  [49]  and  [50]  can  be  equated  to  eliminate  qw" ,  and  yield  a 


differential  equation  for  o  = 


3/7 

1/7  ,, 

/ 

n  d6  * 

+  .*  8  „  do 

- 

72 

dr* 

u  7  dr* 

1/4 

r  2  o  7  3//8 

r* 

Ll  +  6  (  1  -  r*  )  J 

0.0225 


2  IT 


,nr,f 


1/4 


*1/4  1/7 

5  o 

1/7 

Multiplying  [51a]  by  o  ,  and  recognizing  that 

9/7 


[51a] 


2/7 

a  do 


do 


dr*  9  dr* 

allows  [51a]  to  be  written  as: 
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3/7 


9/7 


dr* 


1/4 


1 

d  6  * 

— 

81 

6  * 

dr* 

7 

9 

rj  2  “ 

,  3/8 

1+6 _ (  1  -  f 

*2)  _ 

1 

2  Tt 


N 


R,F 


1/4 


[51b] 


1/4 


S/7 

Eq.  [51b]  is  a  first  order  differential  equation  linear  in  o 
It  is  solved  with  the  aid  of  the  integrating  factor 


J _ d± 

6*  dr* 


dr* 


9/8 


With  the  initial  condition  that  at  r  =  rD  (r*  -  1) ,  6 


0 ,  the 


solution  for 


is: 


3/7 


-  0.413 


<5*\y/5)S/8 


n  1//4  [l  +  B  2(l-n  2)2]  dn. 


J 

I 


(«*  »r,F1/5>1/8 


[52] 


This  is  an  explicit  solution  for  a  as  a  function  of  r*  and  3 

only  (6*  Nr  Is  a  function  of  r*  and  3  only).  The  integral 

in  [52]  must  be  evaluated  numerically  for  all  values  of  6  ,  with 

the  values  of  5  *  N  determined  from  the  solution  of  [ 4 Q]  . 

R,  F 

With  a  evaluated,  the  heat  transfer  from  the  rotating  disk  can  be 
-evaluated  from  either  [49]  or  [50]  .  Eq.  [49]  is  preferable 
as  [50]  involves  differentiation  of  o  .  This  is  considerably 
less  accurate  in  a  numerical  solution  than  raising  °  to  the  l/7th 
power . 

The  solution  of  [52]  yields  values  of  the  boundary  layer 
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ratio  o  =  A  */6  *  greater  than  unity.  0  is  1  at  r*  =  0  ,  on  the 
order  of  5  at  r*  =  0.5,  and  increases  indefinitely  as  r*  decreases 
toward  zero.  This  contradicts  the  implicit  assumption  of  [42] 
where  A  *  is  assumed  to  be  less  than  6  *  so  that  the  l/7th  power 
velocity  profile  applies  over  the  range  of  integration.  Eq.  [  42  ] 
could  be  corrected  by  integrating  from  0  to  6  *  using  the  1/7 th 
power  velocity  profile  and  from  6*  to  A*  using  the  velocity  u 
A  nonlinear  differential  equation  in  o  would  then  result, 
and  a  solution  would  have  to  be  obtained  numerically.  The  extra 
effort  was  not  felt  worthwhile  as  the  solution  is  approximate  and 
based  on  some  unverified  as  sumptions .  Also,  as  o  appears  to 
the  l/7th  power,  large  errors  in  evaluating  °  will  not  affect  the 
heat  transfer  proportionately. 

Defining  a  local  Stanton  number, 


N 


St,  1  = 


w 


1 


^w  ^  A  ^ 


W  '  c 
2*  r  z,  P 


o  “o 


[  49  ]  can  be  written  as: 

0.2 

NSt,l  NR,F 


2  ,  2  1  3/8 

0.0357 

*~S  /  A 

_1  +  6(1-  r*2)  I 

(6*  NRy/5)l/4  a  1/7 


1 5  3  ] 


In  this  analysis,  it  has  been  assumed  that  eTI  =  e 

H  M 

everywhere.  This  implies  that  Npr  =  1.  However,  a  Npr  effect 

can  be  included  if  it  is  assumed  that  the  effect  of  N  on  the 

Pr 
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rotating  disk  heat  transfer  is  the  same  as  on  turbulent  flat  plate 
heat  transfer.  This  is  logical,  as  both  flows  are  essentially  tur¬ 


bulent  boundary  layer  flows.  For  the  Npr  range  of  gases,  the 

0 . 4 

effect  of  Npr  is  felt  as  Np^  .  "Correcting"  [53  ]  for  the 


Np^  effect  yields  the  final  relation  for  heat  transfer: 


0.4  l 
NSt,lNPr  NR,F 


'•2  _  0.03  57  El  +  6  2  (1  -  r*2)2_ 


2  1  3/8 


r*3/4  (6*  NRfF 


l/5)1/ 4  0  1/7 


For  air,  NRr  =  0.7,  and  [54]  becomes: 


NSt,l  Nr,F 


•2  =  0.0411  [l  +  B  2  (1  -  r*2)2  E  3//8 

"  r*V’  <5*  nr,f'/5>1/4  »  1/7 


This  defines  the  local  heat  transfer  from  a  rotating  disk. 

The  right  hand  side  is  a  function  of  r*  and  NR  j^/Nr  R  only. 

1/5 

(6*  Nr  p  )  is  obtained  from  the  solution  of  [40]  ,  and  o  from 
the  solution  of  [52], 

In  Fig .  38 ,  N  N  3  ‘  4  ND  is  presented  as  a  function 

St  Pr 

of  r*  for  various  NR  jy[/NR  R .  The  heat  transfer  is  seen  to  de¬ 
crease  at  first,  and  then  increase  indefinitely  as  the  radius  de¬ 
creases.  The  initial  decrease  is  similar  to  that  for  a  flat  plate; 
as  the  hydrodynamic  boundary  layer  thickens,  the  heat  transfer 
decreases.  However,  as  the  fluid  flows  radially  inward,  the 
fluid  accelerates,  which  thins  the  hydrodynamic  boundary  layer 
and  increases  the  heat  transfer  coefficient.  The  heat  transfer 
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coefficient  increases  indefinitely  near  the  disk  center  as  here  the 

fluid  velocity  becomes  indefinitely  large. 

The  average  heat  transfer  is  obtained  by  integrating  the  local 

heat  transfer  coefficient  over  the  disk  area.  The  average  Stanton 

number  is  thus  related  to  the  local  Stanton  number  by: 

I 


NSt  NPr 


°-4  0.2  (  0.4  0.2 

(NSt,lNPr  NR,F  >  > 


N 


R,F 


[56] 


Eq .  [54  ]  is  used  to  evaluate  the  integrand,  which  is  a  function  of 
r*  and  ^  only.  The  integration  is  performed  numerically. 

At  the  limits  of  integration,  Ngt  ^  becomes  infinite.  At  the 
upper  limit,  the  difficulty  is  circumvented  by  simplifying  [54  ]„ 

For  r*  close  to  unity,  some  of  the  terms  in  [54]  are  small  relative 
to  the  remaining  terms  and  can  be  neglected.  Eq.  [56]  can  then 
be  analytically  integrated  over  the  range  r*  =  0 . 99  to  1 . 00  with  an 
accuracy  of  about  0.1  per  cent.  At  the  lower  limit,  the  effect  of 
the  singularity  is  small  for  ^  =  0,  and  the  integral  can 

be  integrated  numerically.  For  j^/N^  ^  greater  than  zero, 
the  value  of  the  integral  from  an  arbitrary  lower  limit  to  r*  =  1, 
normalized  by  the  value  of  the  integral  over  the  same  range  for 
NR,M//NR,F  =  can  numerically  evaluated  and  plotted  as  a 
function  of  the  lower  limit.  Extrapolation  to  r*  =  0  allows  the 
evaluation  of  the  average  heat  transfer  behavior. 

The  average  heat  transfer  results  obtained  in  this  manner 
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are: 


TABLE  III 

ANALYTICAL  RADIAL  INFLOW  HEAT  TRANSFER 
RESULTS 


m/Nr\f 

0.4  0.2 

NSt  NPr  NR,F 

. .  _ T  0.2 

NSt  Nr,  f 

(NPr  =  0.7) 

0 

0.0S48 

0.1089 

0.1 

0.1018 

0.1169 

0.3 

0.1368 

0.1570 

0.5 

0.1768 

0.2030 

0.6 

0.2000 

0.2298 

0.7 

0.2183 

0.2507 

The  average  heat  transfer  results  are  presented  in  Figs.  30, 
31  and  32  in  comparison  with  the  experimental  results.  The  dis¬ 
crepancies  between  analysis  and  experiment  are  discussed  in 
Section  VII . 

D .  Effect  of  Shroud- to- Disk  Spacing  on  Heat  Transfer 

If  the  disk- to- shroud  spacing  is  finite,  the  boundary  layers 
on  the  disk  and  shroud  will  increase  the  velocity  of  the  fluid 
external  to  the  boundary  layers.  The  boundary  layer  growth  and, 
consequently,  the  heat  transfer  will  be  affected.  The  amount  of 
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the  velocity  increase  is  determined  from  continuity  considerations. 
The  effect  on  the  heat  transfer  can  be  estimated  from  the  increased 
velocity . 

The  velocity  increase  is  evaluated  from  the  continuity  equation. 


The  flow  rate  is  given  by: 


rzo 


W  = 


2  n  r 


u  dz 


J 

O 

=  —  2  tt  r  p 


r  ^  disk  zG 

J  udz  +  U6  (zQ  -  6  disk  -  <5  shroud)  +|  u  dz 

-  °  ^shroud 

Assuming  that  6  ^  ^  =  6  shroud  =  5  '  as  t^ie  boundary  layer  thick¬ 


nesses  are  small,  and  using  the  l/7th  power  velocity  profiles  inside 
of  the  boundary  layers,  the  flow  rate  becomes: 


6  * 

w  =  -  2  n  p  rr0u5  (z0‘-  —  ) 

The  velocity  of  the  fluid  external  to  the  boundary  layer  is  then: 

_ _ W _ 

U  6  2  ,p  r  rQ  (zD*  -  JJ_)  " 

The  velocity  assumed  in  the  analysis  is: 


[  5T  ] 


[58] 


W 


Z  *  =  00 

o 


2  *P  r  ro  V 


The  ratio  of  the  two  velocities  is: 
u  1 


6 


6 


(  1 


7  *  = 


4  zo*' 


[59] 


[60] 


The  effect  of  the  shroud  spacing  is  largest  at  the  large  6  *  which 


97 


occur  at  low  Np  M/Np  y  anc*  at  t'ie  smaH  disk-to-shroud 
spacings . 

The  effect  of  the  velocity  increase  on  heat  transfer  can  be  esti¬ 
mated.  The  heat  transfer  coefficient,  h,  is  proportional  to  Ng^Np  p. 

r  ,  0.8 

From  [  54  J  ,  N„  Nn  is  proportional  to  ND  „ 

bl  I'j  T  I\,l" 


0.8 


Therefore,  h  is 

proportional  to  Nn  „  .  An  increase  in  velocity  due  to  boundary 

layer  growth  can  be  thought  of  as  an  increase  in  velocity  due  to 

increased  flow  rate  (Nn  J) .  The  effect  of  the  shroud-to-disk 

K,  r 

spacing  on  the  heat  transfer  can  then  be  estimated  as: 


0.8 


N 


h) 


St 


N 


r7  *  _  00 


st) 


u6) 


7  *  =  oo 


[61] 


where  u  &  /ufi  )  Zq* 


is  evaluated  from  [  60  ] 


For  z  *=0.1  and  N  =10  (the  minimum  values  to  be  studied 

°  R,  F 

experimentally),  and  for  N„  ../N  =0,  the  increase  in  heat  trans- 

f  JV1  K  j  r 

fer  is  about  3  per  cent  at  most.  Because  5*  decreases  at  higher 
values  of  Nr  p'  effect  will  decrease  as  NR;M/Np  p  in¬ 

creases.  This  leads  to  the  important  conclusion  that  zQ*  may  not 
be  an  important  parameter  in  the  flow  and  heat  transfer  studies. 
Elimination  of  zQ*  as  a  parameter  would  considerably  reduce  the 
length  of  the  experimental  investigation.  This  conclusion  must  be 
verified  experimentally. 
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E .  Conclusions 


The  conclusions  that  can  be  drawn  from  the  results  of  the  ana¬ 
lytical  investigation  are: 

1.  The  problem  of  heat  transfer  from  a  shrouded  rotating  disk 
has  been  attacked  using  integral  techniques.  The  relation 

-  —  between  the  nondimensional  parameters  has  been  predicted. 

There  is  considerable  disagreement  between  the  analytical 
and  experimental  fluid  dynamics  and  heat  transfer  behavior. 

2.  The  predicted  effect  of  the  shroud-to-disk  spacing,  zQ* ,  on 
the  heat  transfer  is  small.  The  predicted  heat  transfer  for 
zq*  =  0.1  is  at  most  only  3  per  cent  higher  than  that  for 
zo*  =  "  •  The  lack  of  an  effect  of  the  disk-to- shroud 
spacing  will  have  to  be  verified  experimentally. 

3.  The  present  analytical  results  will  apply  to  a  rotating  disk 

with  closely  spaced  vanes.  If  it  is  assumed  that  the  vanes 
prevent  any  relative  tangential  movement  (v  =  no  every¬ 
where),  the  same  solution  as  for  NR  ,,/N  =  0  is  ob- 

tained.  However,  on  a  vaned  rotating  disk  there  would 
probably  be  a  secondary  flow,  although  not  of  the  same 
magnitude  as  found  on  the  unvaned  disk.  The  effect  of 
vanes  on  the  heat  transfer  behavior  will  have  to  be  deter¬ 
mined  experimentally. 
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FIGURE  3  6 


PREDICTED  FLOW  PATHS  WITH  RESPECT 
TO  THE  DISK 


FIGURE  3  7 

ANALYTICAL  HYDRODYNAMIC  SOLUTION 

6*  Nd  °-2  AS  A  FUNC'TION  OF  r* 

K,  t 


FIGURE  38 

ANALYTICAL  HEAT  TRANSFER  SOLUTION 
NSt  lNPr°*4NR  F°'2  AS  A  FUNCTION  OF  r* 
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VII  DISCUSSION  OF  RESULTS 


The  analytical  and  experimental  heat  transfer  results  differ  signi¬ 
ficantly  in  the  magnitudes  of  the  heat  transfer  coefficients  and  of  the 
Reynolds  number  exponents.  The  discrepancies  are  due  to  some  over¬ 
simplifications  made  in  the  analysis.  The  assumptions  will  be  exam¬ 
ined  in  light  of  the  experimental  results. 

The  assumption  of  axial  symmetry  is  experimentally  verified  for 

N  a/Nd  p  greater  than  0 . 3  .  The  assumptions  concerning  the 

details  of  the  boundary  layer  flow  are  quite  reasonable.  The  boundary 
layer  flow  is  wall  dominated;  there  are  no  wakes  or  jets  to  influence 
the  boundary  layer.  This  implies  that  the  "law  of  the  wall"  (l/7th 
power  velocity  and  temperature  profiles)  is  valid. 

The  most  questionable  simplifications  are  those  concerning  the 
influence  of  the  disk  and  shroud  boundary  layers  on  the  flow  pattern. 
The  assumed  influence  is  an  increase  in  the  one-dimensional  fluid 
velocity  external  to  the  boundary  layers  due  to  the  reduced  veloci¬ 
ties  in  the  boundary  layers.  The  actual  effect  is  that  the  shear 
stresses  on  the  disk  and  shroud  reduce  the  fluid  angular  momentum 
near  the  two  surfaces.  The  resulting  unbalance  in  centrifugal  force 
field  across  the  flow  passage  forces  the  fluid  to  flow  radially  out¬ 
ward  in  the  passage  center,  as  shown  in  Fig.  20.  This  influence 
of  the  surfaces  is  similar  to  that  on  the  flow  between  rotating  con¬ 
centric  cylinders,  in  which  a  system  of  counterrotating  vortices 
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unexpectedly  is  formed,  (24). 

The  flow  pattern  is  thus  governed  by  the  presence  of  the  disk  and 
shroud.  The  flow  can  not  be  considered  as  an  external  flow  with 
boundary  layers  on  the  surfaces.  While  the  boundary  layer  profile 
is  quite  probably  of  the  l/7th  power  form,  the  radial  velocity  in  the 


nr,f  w 

i 

boundary  layer  does  not  approach 

1  iv  rQ 

r* 

as  assumed 

(Eq.  [29]  ).  This  is  seen  in  Figs.  14  to  18,  where  the  external  one¬ 
dimensional  velocity  does  not  characterize  the  velocity  profiles. 

This  breakdown  of  the  concept  of  an  external  flow  with  boundary 
layers  is  reinforced  by  the  comparison  between  the  experimental  and 
the  predicted  heat  transfer.  The  heat  transfer  exponent  is  directly 
related  to  the  velocity  profile  power  in  a  boundary  layer  analysis. 

A  1/n  power  profile  yields  a  2/(n  +  3)  power  on  Reynolds  number; 

0.2 

a  l/7th  power  yields  Np  p  .  The  velocity  profile  power  corre¬ 
sponding  to  the  experimentally  determined  0.55  power  on  Np  p  would 
be  1/0.638.  This  order  of  magnitude  difference  in  exponent  and 
the  reversed  curvature  of  the  profile  are  physically  unrealistic. 


l/7th  power  profile  1/0.638  power  profile 
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The  heat  transfer  coefficients  are  based  on  the  measured  tem¬ 
perature  difference  between  the  disk  and  the  air  stream  inlet.  At 

hiqh  N„  ,,  and  Nn  the  disk  is  rotating  rapidly  and  the  rela- 
R f  Rul  K/  r 

tive  velocities  are  quite  high.  As  a  result,  the  difference  between 

o,  /  r 

the  adiabatic  wall  temperature  and  the  d^sk  temperature  is  signifi¬ 
cant  for  some  test  conditions.  The  difference  between  the  two 
temperatures  is  given  by: 


adiabatic  wall  air 


-  t  •  )  =  Nt 

sir  '  \ 


2  gc  C 


1/3 

where  Npr  =0,9  is  the  recovery  factor  for  turbulent  radial  out¬ 
flow  over  a  rotating  disk  (2  3).  As  the  fluid  relative  velocity  VR 
is  a  function  of  radius,  the  temperature  difference  as  a  function  of 
radius  must  be  evaluated  and  an  area  average  value  obtained. 

The  maximum  computed  average  difference  between 

t  ,,  and  t  .  is  1.8  F,  increasing  the  heat  transfer 

adiabatic  wall  air 

coefficient  by  about  4  per  cent  over  the  reported  values.  However, 

2  i2 

as  the  temperature  difference  decreases  both  as  (NR  M)  and(NR,R) 
decrease,  the  effect  on  the  heat  transfer  coefficients  decreases 
rapidly.  For  about  90  per  cent  of  the  heat  transfer  data,  the  effect  on 
the  heat  transfer  coefficient  is  less  than  1  per  cent.  Therefore,  the 
data  was  not  corrected  for  this  effect. 

The  fluid  dynamics  investigation  was  conducted  with  the  disk 
and  air  stream  at  the  same  temperature,  while  in  the  heat  transfer 
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studies  a  temperature  difference  existed.  The  flow  patterns  in  these 
two  studies  were  possibly  different  due  to  the  free  convection  effect 
of  the  centrifugal  force  field  on  the  density  differences  present  in 
the  heat  transfer  study. 

The  parameter  characterizing  the  free  convection  is  the  Grashof 
number: 


For  the  disk  tests,  the  acceleration,  g,  is  the  centrifugal  accelera- 
2 

tlon  ro  “  '  Using  the  ideal  gas  relation  for  the  coefficient  of 

thermal  expansion,  8  =  1/T,  the  Grashof  number  can  be  written  as: 


over  the  rotating  disk  test  range . 

The  parameter  characterizing  the  through  flow  in  the  same 
manner  as  does  the  length  Reynolds  number  for  flow  over  a  flat 
plate  is: 


2  IT 


The  flow  produced  by  the  free  convection  effect  is  in  the 
same  direction  as  the  through  flow.  As  the  centrifugal  force  acts 


radially  outward,  the  less  dense  (higher  temperature)  air  next  to  the 
disk  is  forced  radially  inward.  Eckert  and  Diaguila  have  established 


the  following  criteria  for  the  combined  effects  of  free  and  forced  con¬ 
vection  on  flat  plate  heat  transfer  (5): 

For  a  deviation  in  the  heat  transfer  of  less  than  10  per  cent  from 
pure  free  convection: 


NGr 

0.007  (NRe)2-5 


> 


1 


or 


1300 

OTT 


I  Nr,m\2>  1 

\NR,F  / 


[63] 


For  a  deviation  in  the  heat  transfer  of  less  than  10  per  cent 


from  pure  forced  convection: 

0.00  16  (NRe)2'5 


>  1 


N, 


Gr 


or 


0.00018  NR/F°-5  >  1 


[6M 


Evaluating  these  criteria  over  the  heat  transfer  test  range  yields 
the  following  flow  divisions: 


NR,M//Nr,F 


0 . 31 

0 . 66 


0.99 


Type  of  Flow 
Forced  convection 
Forced  and  free  convection 
Free  convection 
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At  high  NR,M/NR  F'  the  flow  pattern  would  be  expected  to  be 
dominated  by  the  free  convection  current.  However,  an  attempt  to 
predict  the  disk  heat  transfer  using  flat  plate  convection  correlations 
was  unsuccessful.  The  predicted  free  convection  coefficients  are' 
one- third  to  one-half  those  measured,  and  differ  considerably  in 
their  dependence  on  NR  ^  and  NR  p.  This  lack  of  correspondence 
at  high  NR  jyj/NR  p  may  imply  that  the  flat  plate  free  convection 
criteria  are  not  applicable  to  the  rotating  disk  flow. 

The  free  convection  effect  on  the  present  heat  transfer  results 
is  probably  small.  The  heat  transfer  results  of  Fig.  31  show  essen¬ 
tially  the  same  dependence  of  Ng^  on  NR  p  over  the  entire 
Nr.m/^r  p  range.  This  would  not  be  expected  if  regimes  of 
forced  convection,  forced  and  free  convection,  and  free  convection 
occurred  as  NR  m/Nr  p  increased.  At  high  NR  m/Nr  p,  the 
through  flow  Reynolds  number  NR  p  would  then  not  be  significant, 
and  the  heat  transfer  data  would  then  be  correlated  by  the  Grashof 
number,  or  NR  M.  The  present  correlation  of  the  heat  transfer 
data  over  the  test  range  of  NR  j^/N^  p  strengthens  the  acceptance 
of  the  present  radial  inflow  heat  transfer  results  as  representative 
of  rotating  disk  heat  transfer. 


106 


VIII  CONCLUSIONS 


The  following  conclusions  can  be  drawn  from  the  results  of  the 
present-  investigation: 

1.  A  test  facility  has  been  constructed  that  will  allow  deter¬ 
mination  of  the  average  heat  transfer  from  a  rotating  disk. 

The  facility  is  adaptable  to  studies  of  various  flow  con¬ 
figurations  . 

2  .  The  flow  pattern  for  a  rotating  disk  with  a  single  radial 

inflow  has  been  experimentally  determined.  The  presence 
of  the  disk  and  shroud  have  been  shown  to  govern  the 
fluid  behavior. 

3.  The  heat  transfer  from  a  rotating  disk  with  a  single  radial 
inflow  has  been  experimentally  determined.  At  a  disk-to- 
shroud  clearance  ratio  of  zQ*  =  0.113  and  over  a  Np  p 
range  of  1 . 5  x  10  5  to  6 . 6  x  10  5  and  a  NR  ^j/Nr  p  range 
of  0.1  to  1.7,  the  heat  transfer  is  correlated  by 

NStNR,F°'55'  23’81  +  4'°68  (NR.m/Nr,F)  1111 

4.  The  analytical  investigation  did  not  adequately  predict 
the  fluid  dynamics  and  heat  transfer  behavior.  The 
assumptions  regarding  the  interaction  between  the  disk 
and  shroud  boundary  layers  and  the  main  flow  were  over¬ 
simplified  . 
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IX  RECOMMENDATIONS  FOR  FURTHER  WORK 


The  following  recommendations  for  further  work  are  made: 

1.  The  test  facility  should  be  improved  to  eliminate  the  flow 
leakage  past  the  rotating  screens.  The  possibility  of 
modifying  the  facility  to  permit  the  determination  of  local 
heat  transfer  coefficients  should  be  considered. 

2.  The  effect  of  free  convection  on  the  rotating  disk  heat 
transfer  should  be  determined.  This  would  necessitate 
obtaining  velocity  profiles  during  heat  transfer  tests. 

3.  An  investigation  of  the  details  of  the  boundary  layer  flow 
on  the  disk  should  be  performed.  This  study,  coupled 
with  a  more  realistic  analytical  approach,  should  yield 
more  insight  into  the  fluid  behavior. 

4.  Other  rotating  disk  flow  configurations  should  be  studied 
in  an  approach  to  the  modeling  of  the  veil-cooled  radial 
flow  gas  turbine.  A  logical  sequence  is: 

a)  The  rotating  disk  with  a  veil  flow  and  a  main  flow. 

b)  A  vaned  rotating  disk  with  a  single  radial  inflow. 

c)  A  vaned  rotating  disk  with  a  veil  flow  and  a  main 
flow . 
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X  APPENDICES 


Appendix  A  Experimental  Uncertainties 

The  following  uncertainties  are  considered  representative  for 
the  rotating  disk  test  range. 

1 .  Dimensions 

At  20:1  odds,  the  uncertainty  intervals  are  estimated 
to  be: 


ro 

7.979" 

+ 

0.001" 

+ 

0.01% 

z 

0.900" 

+ 

0.010" 

+ 

1.1% 

o 

Disk  area 

+ 

0.02  % 

2 .  Fluid  Properties 

At  20:1  odds,  the  uncertainty  intervals  for  the  fluid  pro¬ 
perties  are  estimated  to  be: 


T 

540  R  - 

0. 

1  R 

O 

+  1 

.02% 

P 

30"Hg  - 

0. 

0  5"Hg 

o 

+  1 

.2% 

p  ±  0.2% 

Cp  ±  0.2% 

V  -  0.2% 

k  -  0.3% 

3  .  Air  Flow  Rate  ,  W 

The  air  flow  metering  section  was  constructed  according 
to  the  specifications  of  the  A.S.M.E.  Power  Test  Code  (1). 
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The  metering  section  was  calibrated  against  a  standard  flow 
metering  tank.  The  metering  section  orifice  coefficient  was 
computed  using  the  flow  rates  determined  from  the  flow  meter 
tank  and  the  pertinent  temperatures  and  pressures  measured 
at  the  orifice.  The  calibrated  coefficient  was  then  compared 
to  that  tabulated  in  (1) . 

At  20:1  odds,  the  flow  rate  uncertainty  interval  is  esti¬ 
mated  to  be: 

W  ±1.6% 

4.  Disk  Angular  Velocity,  a> 

The  disk  angular  velocity  was  measured  with  a  General 
Radio  Strobotac  calibrated  to  -  1%.  The  variation  in  disk 

speed  during  a  test  is  estimated  to  be  t  1%.  At  20:1  odds, 
the  disk  angular  velocity  uncertainty  interval  is  estimated 
to  be: 

u  t  1 . 5  % 

5 .  Heat  Flux  Through  the  Disk,  q" 

The  uncertainty  in  the  heat  flux  passing  through  the  disk 
is  estimated  from  computations  of  the  disk  thermal  resistance. 
The  thermal  resistance  is  defined  as  the  temperature  difference 
between  the  heat  meter  and  test  surface  thermocouples  divided 
by  the  heat  flux.  The  thermal  resistance  was  continuously 
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evaluated  during  calibration  and  testing.  The  variation  in 
thermal  resistance  is  assumed  to  be  due  solely  to  heat  flux 
measurement  uncertainties. 

At  20:1  odds,  the  heat  flux  uncertainty  interval  is  esti¬ 
mated  to  be: 

q"  ±  2.0% 


6 •  Temperature  Difference  Between  Disk  Test  Surface  Anri 
Air  Stream 

The  copper-constantan  thermocouples  were  calibrated 
against  primary  standard  thermometers  in  a  controlled  tempera¬ 
ture  oil  bath.  At  20:1  odds,  the  air  temperature  uncertainty  is 
estimated  to  be: 

W  -  0.1  F 

The  uncertainty  in  a  single  disk  surface  temperature 
measurement  due  to  extraneous  EMF's  in  the  garter  spring 
pickup  and  to  uncertainties  in  the  calibration  is  ir  0.3  F. 

The  test  surface  temperature  variation  is  on  the  order  of 

-  1.0  F.  The  estimated  test  surface  temperature  uncertainty 
is  then 


disk 


t  1 .3  F 


At  20:1  odds,  the  estimated  uncertainty  of  the  tempera¬ 
ture  difference  between  the  disk  and  the  air  stream  is 
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estimated  to  be: 


t  -  t 
disk  air 


50  F  1 


=  ±2,8% 


7 .  Nondimensional  Parameters 

At  20:1  odds,  the  estimated  uncertainties  in  the  non- 
dimensional  parameters  are: 


NNu 

+ 

3.4% 

nr,f 

± 

1.7% 

nr,m 

± 

1.5% 

Nst 

± 

3.9% 

NStNR,F°'55 

± 

4.0% 
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Appendix  B  Fluid  Velocity  Mea surement  Errors 

The  velocity  and  flow  angle  measurements  are  subject  to  a 
number  of  possible  errors.  The  greatest  source  of  error  in  the 
velocity  determination  stems  from  the  method  of  measuring  the 
dynamic  pressure.  The  velocity  probe  is  equipped  to  sense  total 
pressure  only;  the  static  pressure  must  be  measured  separately. 
With  the  probe  mounted  in  the  flow  passage,  the  wall  static  pres¬ 
sure  can  not  be  measured  at  the  same  location  as  the  total  pres¬ 
sure.  Therefore,  to  obtain  a  dynamic  pressure  measurement,  the 
following  procedure  was  used. 

The  total  pressure  tap  was  connected  to  one  leg  of  a  differen¬ 
tial  manometer.  A  wall  pressure  tap  at  the  same  radius  as  the 
probe,  but  at  a  different  circumferential  position,  was  connected 
to  the  other  leg  of  the  manometer.  A  velocity  traverse  across  the 
flow  passage  was  then  made,  and  the  pressure  difference  recorded. 
At  the  same  test  conditions,  the  difference  in  static  pressure  at 
the  wall  between  the  static  and  total  pressure  measuring  locations 
was  measured.  The  dynamic  pressure  was  then  corrected  for  this 
difference  so  that  the  total  and  static  pressure  were  determined  at 
the  same  location.  The  correction  varied  between  1  and  30  per  cent 
of  the  dynamic  pressure.  The  resulting  uncertainty  in  the  velocity 
is  estimated  to  be  -  3  per  cent. 

The  uncertainty  in  the  dynamic  pressure  is  further  increased 
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in  that  the  static  pressure  could  not  be  measured  at  the  same  point 
in  the  flow  passage  as  the  dynamic  pressure.  As  seen  in  Figs.  14  to 
18,  the  velocity  variation  across  the  passage  is  large,  and  a  static 
pressure  variation  would  also  be  expected.  No  estimation  of  the 
static  pressure  variation  was  made. 

The  large  velocity  variations  across  the  flow  passage  cause 
a  curving  of  the  streamlines  in  the  vicinity  of  the  probe.  The  higher 
velocities  shift  toward  the  probe  tip;  the  probe  measures  a  velocity 
higher  than  actually  exists  in  the  absence  of  the  probe.  The  fluid 
displacement  is  on  the  order  of  0 .02" ,  but  is  not  known  accurately 
enough  to  correct  the  traverse  results. 

The  probe  was  calibrated  in  a  wind  tunnel  for  flow  angle  measure 
ments  to-  -i  0.2  5°.  The  uncertainty  in  reading  the  protractor  is 
±  0.25°.  The  uncertainty  interval  for  flow  angle  measurements  due 

to  instrument  uncertainties  is  t  0.5°. 

The  flow  angle  measurements  are  sensitive  to  velocity  gradients 
normal  to  the  probe  axis.  At  high  j^j/N^  F'  ^ow  *s  essen_ 
tially  in  the  tangential  direction,  and  the  velocity  gradient  in  the 
radial  direction  is  large,  as  shown  in  the  following  figure.  Conse¬ 
quently,  the  two  directional  taps,  which  are  separated  by  the  probe 
diameter,  may  be  in  regions  of  different  total  pressure.  At  a  null 
pressure  indication,  the  probe  is  then  orientated  so  that  the  tap 
in  the  low  total  pressure  region  measures  a  higher  component 
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of  the  dynamic  pressure  than  does  the  other  tap.  The  probe  is  thus 
rotated  too  far  in  the  direction  of  the  higher  velocities.  This  flow 
angle  error  is  estimated  to  be  within  2°  for  the  rotating  disk  tests. 


Flow  Direction 

- ► 


The  radial  velocity  uncertainty  is  largest  at  high  NR  Wf- 
The  radial  velocity  is  obtained  from  the  product  of  the  resultant 
velocity  and  the  sine  of  the  flow  angle.  For  NR  ,m/nr,f  greater 
than  0.6,  the  flow  angle  is  between  80°  and  100°.  An  uncertainty 
of  -  0.5°  results  in  a  radial  velocity  uncertainty  on  the  order 
of  t  25  per  cent. 
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Appendix  C  Data  Reduction  Procedure 


A  Burroughs  220  digital  computer  was  used  to  reduce  the  data 
with  the  exception  of  the  mean  heat  transfer  data.  Separate  programs 
were  written  for  the  computation  of  Np  p  and  Np  j^/Np  p,  and  for 
the  determination  of  the  velocity  profiles.  These  programs  are 
written  in  the  Balgol  language. 


1 .  Evaluation  of  Np  p  and  Nr  jyp/Npp 

The  following  table  gives  the  translation  between  the 
present  nomenclature  and  that  used  in  the  computer  program. 


ATM 

D,  O,  or  OR 
I 

DM 

DO 

DPO 

KO 

MU 

MV 

R 

P 

T 

YO 


denotes  atmospheric  conditions 
denotes  orifice  conditions 
denotes  inlet  conditions 
density  correction  factor  for  humidity 
orifice  diameter 

pressure  drop  across  the  orifice 

orifice  coefficient 

dynamic  viscosity 

temperature  in  millivolts 

density 

pressure 

temperature  in  F 
orifice  expansion  factor 
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MFLO 

W,  mass  flow  rate 

NRF 

NR,F 

NRM 

nr,m 

NRM/NRF 

nr,m/nr,f 

RPM 

u>  ,  disk  angular  velocity 

RUN 

Run  number 

ZO 

zo 

Empirically  determined  relations  for  KO,  MU,  T,  and  YO 


were  employed.  The  machine  program  follows. 


COMMENT  DISK  FLOW  DATA  REDUCTIONS 
INTEGER  RPM,  RUNS 
WRITE  ( SSFMT ) S 

FORMAT  FMTIB5»*MFLO*»B5  ,*RPM*  »B3  »*RHO  ! N*»B5  **NRF*»B10»*NRM**B7  » 
♦NRM/NRF*  *B3 ,*RUN*.W4 ) $ 

START  o  .READ  (SSDATA)S 

INPUT  DATA ! PATM  »  PD  #P  I *DPO»MVO*MVI » DM » RPM* DO* RUN  )  S 
TOR  *  ( SORT I MVO  +  9.081)  -  2.899 ) /0. 003567$ 

RO  *  1 o3207 ( PATM  +  PD)/ltOR  +  459.7)$ 

MUO  *  0.0399  +  0c0000605tOR$ 

EITHER  IF  DO  GEO  3.0$  GO  BIGS 

OTHERWISES  GO  LlTLS 

BIG..  COMMENT  3.351  ORIFICES 

NRD  *  66201  SORT  C  DPO.RO) ) /MUOS 
LNRD  *  LOG  I NRD)  $ 

KO  *  0*757254  -  0 . 0 1 75707LNRD  +  0 . 00064 5 328 LNRD. LNRD $ 

YO  «=  1*0  -  0.0234DPO/IPATM  +  PD»S 
MFLO  -  4033K0*  YO.DM. SORT  <  DPO.RO) S 
GO  OUTS 

LI TL. .COMMENT  1.919  ORIFICES 

NRD  »  20341  SORT (DPO.RO) ) /MUOS 
LNRD  *  LOG i NRD ) $ 

KO  «  0.692822  -  0.0152481LNRD  +  0 . 000622632LNRD.LNRDS 
YO  =  1.0  -  0.0218DPO/M PATM  +  PD) $ 

MFLO  a  1321KO.YO.DM.SQRTI DPO.RO ) $ 

OUT..ZO  «  0.07580$ 

TI  «  li  SORT  (  MV  I  +  9.081)  -  2 . 899  ) /0* 003567S 
MU  I  *»  0.0399  +  0.0000605TIS 
NRF  *  MFLO/ <  ZO. MUD  $ 

RI  «  1 o  3207 ( PATM  +  PI)/(TI  +  459.7)$ 

NRM  »  1 6 6 • 7 R PM .RI/MUIS 

11.8 


WRITE  <**ANS»HDG)$ 

OUTPUT  ANSI MFLO*  RPM  *R I ,NRF • NRM ♦NRM/NRF »RUN » $ 
FORMAT  HDG<X11.1.I7,X9.5»F13.5»F13.5*X9.4»16.W4}S 
GO  STARTS 
FINISHS 
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2  .  Evaluation  of  Velocity  Profiles 


The  following  table  gives  the  translation  between  the  present 
nomenclature  and  that  Used  in  the  computer  program: 


A 

ANG 

AO 

DPV 

DPVC 

MW 

PATM 

PV 

ROV 

TV 

V 

VN 

VT 


measured  flow  angle 

flow  angle  relative  to  the  radial  direction 

protractor  angle  for  radial  flow 

dynamic  pressure  as  measured 

static  pressure  correction 

temperature  in  millivolts 

atmospheric  pressure 

static  pressure 

density 

temperature  in  F 
resultant  velocity 
radial  velocity  component 
tangential  velocity  component 


The  machine  program  follows. 
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COMMENT  VELOCITY  PROFILES 
INTEGER  I, N, RUNS 
ARRAY  DPVI 20 ) »  A«20>$ 

START ..READ (SSOATAjS 

INPUT  DATAtRUN»PATM.PV#MVV,AO»N,DPVt, 

FOR  I  *  H«i.NI$DPVU  I  »  FOR  I  »  (  1  *  1  *N  J  SA  <  I  )  i  S 
TV  «  (  SORT ( MW  +  9,081  »  -  2. 899  ) /O. 003567$ 

ROV  -  1 • 3207 ( PATM  +  PVI/ITV  +  459o7»$ 

WRITE  (SSANS1»HDG1 JS 

OUTPUT  ANSI (RUN.ROVIS 

FORMAT  HDG1 IB5«*RUN*» 1 4  s B5  **RH0  *«  ,X7 • 5 ,B6 » *VEL* »B5 • *VTAN* ,B5 « 

♦VRAD  *  »B4 »* ANGLE*  »  W4  )  $ 

FOR  I  *  (1,1. N»$ 

BEGIN 

V  «  SORT!  334,88  C0PV«  i'>  +  DPVO/ROVKS 
ANG  *  All)  -  AO  -  3.25$ 

VN  *  V,COS<!  0,0 1745 3 3 ANG)  $ 

VT  «  V «S I N ( 0«0174533ANG!$ 

WRITE  «S$ANS2»H0G2>$ 

OUTPUT  ANS2 ( V  *  VT,  VN »  ANG  S  $ 

FORMAT  HOG 2  I B33 »X6 . 2 , B3 » X6 . 2 » B3 * X 7 . 2 » B3 , X6 . 2 » W4 ) $ 

ENDS 

GO  STARTS 


FINISHS 


3  •  Evaluation  of  Average  Heat  Transfer 


The  average  Stanton  number  was  obtained  from. hand  calcu¬ 
lations.  The  heat  meter  temperatures  were  measured  by  the 
heat  meter  thermocouples.  The  heat  meter  calibration  curves 
(Appendix  E)  were  entered  to  obtain  the  heat  meter  calibration 

factor  in  BTU/hr  ft  _  The  product  of  the  heat  meter  output 
mv 

in  millivolts  and  the  calibration  factor  determined  the  heat  flux 

passing  through  the  heat  meter.  These  heat  fluxes  were  plotted 
2 

as  a  function  of  r*  (proportional  to  disk  area)  and  the  area 
average  heat  flux  was  obtained  using  a  planimeter. 

The  average  disk  temperature  was  obtained  as  the  numeri¬ 
cal  average  of  the  five  test  surface  thermocouple  temperatures. 
This  method  was  sufficiently  accurate  as  the  test  surface  thermo¬ 
couples  are  located  in  the  center  of  equal  disk  area  sections  and 
the  temperature  variation  across  the  surface  was  within  t  1  F. 
The  air  inlet  temperature  was  taken  from  the  averaged  readings 
of  the  four  thermocouples  equally  spaced  around  the  circumfer¬ 
ence  of  the  air  inlet  passage. 

The  average  heat  flux  was  divided  by  the  temperature  dif¬ 
ference  between  the  disk  test  surface  and  the  inlet  air  to  deter¬ 
mine  the  average  heat  transfer  coefficient.  The  average  Stanton 
number  is  defined  as: 
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N 


Q"  /  (  tdisk  -  tair  in  ) 


St 


( 


W 


■) 


2  11  ro  zo 


Appendix  D  Local  Heat  Flux  Data  Reduction 


The  local  Stanton  numbers  are  obtained  from  the  local  values  of 
test  surface  heat  flux.  The  test  surface  heat  flux  can  be  related  to 
the  heat  meter  heat  flux  and  the  test  surface  temperature  distribution 
by  the  solution  for  heat  conduction  through  the  disk.  As  this  solu¬ 
tion  for  the  rotating  disk  boundary  conditions  is  not  available  in  the 
literature,  it  will  be  presented  here. 

The  coordinate  system  and  notation  is: 


For  steady  axisymmetric  heat  flow  with  isotropic  thermal  con¬ 
ductivity,  the  conduction  equation  is 


3  2  t 


1 


3  t 


_L  =  o 


[Dl] 


3r‘  r  3  r  3 

The  boundary  conditions  corresponding  to  the  rotating  disk  test 
conditions  are: 


t  (r,  o)  =  tQ  (rS-)  =  tQ  (  r*  ) 


[D2] 


(  t  (r*)  is  the  experimentally  determined  test  surface  temperature 
distribution),  and 

at  (r ,  b)  _  -  q"b  (“r^)  _  -q"b  (r*) 

3z  k  k 


[D3] 
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(<?"b  (r*)  is  the  experimentally  determined  heat  meter  heat  flux  dis¬ 
tribution)  ,  and 


— 5  t  (roVz)  =  0 

9  r 

(there  is  no  heat  transferred  through  the  outer  edge  of  the  disk)  ■ 
Assuming  a  product  solution  for  t(r,z), 
t  (r ,  z)  =  R(r)  •  Z(z) 

the  variables  R(r)  and  Z(z)  may  be  separated,  and  two  ordinary 
differential  equations  result. 

R"  +  _J_  JV_  _  Z"  =  _  B2 

R  r  R  Z 

2 

where  S  is  a  constant.  Or: 

R"  +  — ■ -  R'  +  6^  R  =0 

r 

Z"  -  B  2  z  =  o 


[D5] 

[D6] 

[D7] 


The  solution  to  [  D  6  ]  for  a  disk  without  a  hole  in  the  center  is  a 
Bessel  function  of  the  First  Kind  and  zero  order  [  (15),  p.  3S  ]  : 

R(r)  =  JQ  (6  r)  [D8] 

The  boundary  condition  at  r  =  rQ  [  D  4  ]  is  applied  to  obtain  the 
eigenvalues.  The  Bessel  function  recurrsion  relation  is  employed 
[  (15) ,  p .  31  ] : 


_d_  Jo  (x) 
dx 


-  Jx  (x) 


[D9] 


Therefore 


_dR  I  =  d  Jq  (r  b  ) 
dr  'ro  d  (r  B) 

ro 

The  relation  [  D  .10  ]  is  satisfied  if 


-  B  Jj  (rQ  B  )  =  0 

(rQ  6)  takes  on  the  values  of 


[DIO] 
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the  zeros  of  (x) .  Denoting  the  zeros  of  (x)  as  a  m  >  the 
allowable  .values  of  6  are: 


CX  rn 

6  =  -  m  =  0,  1,  2,  .  .  .  [Dll] 

ro 

The  solution  [  D  8  ]  is  then: 

R(r)  =  I0  («mf[  )  =  J0  (  am  r*  )  [D12] 

Eq.  [D  7  ]  can  be  directly  integrated: 

Z(z)  =  Ame  ^  r°  +  Bm  e  “m  ro  [D13] 

where  Am  and  Bm  are  constants  of  integration.  The  temperature 
distribution  in  the  disk  is  thus  given  by: 

CO 

t  (r,  z)  =  I  J  (  c»m  r*) 
zn=o 

[DlU] 

This  is  the  general  solution  for  the  axisymmetric  temperature 

distribution  in  a  circular  disk  with  an  insulated  perimeter.  The 

constants,  A„  and  B  ,  are  evaluated  by  applying  the  boundary 
m  m 

conditions  at  the  upper  and  lower  surfaces.  Applying  the  test 
(lower)  surface  boundary  condition  [  D  2  ]: 


m  rr 


-  “ 


m  r. 


A^  e 
m 


+  Bm  e 


oo 

t(r,o)  =  tQ  (r*)  =  l  JQ  (  “mr*)(Am  +  Bm) 


[D15] 


m=o 

To  evaluate  each  Am  and  Bm,  both  sides  of  [D  15]  are  multiplied 
by  r*  J  (  an  r*)  dr*  and  integrated  from  r*  =  0  to  r*  =  1. 

With  n  as  the  dummy  variable  of  integration,  [D  15]  becomes: 
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This  is  one  relation  for  Am  and  B  .  The  other  relation  is 
similarly  obtained  by  applying  the  heat  meter  (upper)  surface 
boundary  condition  [  D  3  ]  : 


Multiplying  both  sides  of  [  D  19  )  by  r*  JD  (  an  r*)  dr*,  inte¬ 
grating  from  r*  =  0  to  r*  =  1,  and  using  the  orthogonality 
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relation  [  D.  17]  ,  yields: 

i 


n  q"b  (  n  )  JD  n  )  d 

k 


L' 


r  b  b~1 

j  “m  r  “m  i- 

—  — 

[_Am  e  °  -  Bm  e 

Jo  <  ) 

[D20] 


The  two  relations  [  D  18]  and  [  D  20]  can  be  solved  for  the  con¬ 
stants  Am  and  separately: 


m 


Am  =  Le 


am^— 

o 


i  n  l0  (n  )  JQ  (  a  n  )  d  n 


Jo  ^  a  m  )J 


cosh  (am  b_  ) 


r>  q"k  (  n  )  JD  (  am  ri  )  d  n 


a  k 
m  ° 


Bm=  Le 


ro 


4  n  t0  ( n  )  Jn  (  am  n  )  d  n 


Jo  (  an  >1  2  cosh  (%  -k) 

-  I  f-y. 


n  (  n  )  J  (  a  n  )  d 


am  k 


o  m 


t  D21  ] 


The  heat  flux  at  the  test  surface  can  now  be  evaluated.  The  test 
surface  heat  flux  is  obtained  from  the  temperature  distribution 
by: 


q"o  <r*) 


k 


5t  (r,o) 
8  z 


[  D22  ] 
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or 


q"0(r*)  -  2  J 


Jo  (  «„  r*  ) 


m=o 


cosh(  am_b_) 

ri0  (am ) 

r 

ro 

L_ 

q"h  (n  ) 


%k  sinh  (an  _b.)  tQ  (n  )  n  JQ(  n  )  dn 
ro  ro  _J 


[D23] 


Recognizing  that  for  m  =  0,  a  =  0  and  JQ  (0)  —  1,  [  D  2  3]becomes: 


I 

q"  (r*)  =  2  Jn  q"  (n  )  d  n  +  2  y 

°  J  b  L 

o  m=o 


Jo  (“mr*  ) 


cosh  (  am^)  [lo  (  am  U 


r. 


J 


q",  ( n  )  -  k  sinh  (  amjj)  t  (  n) 

D  r  r  ^ 


"  JG  (  am  n  )  d  n 


[D2U] 


As  n  ,  the  dummy  variable  of  integration,  is  written  for  r*,  the 
first  term  can  be  written  as: 


d 


(r*) 


2 


which  is  the  area  average  heat  flux  passing  through  the  disk.  The 
remaining  terms,  then,  average  to  zero.  These  two  requirements 
provide  a  check  on  the  numerical  solution. 
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The  Burroughs  220  digital  computer  was  employed  to  solve  [  D  24] 
for  given  heat  flux  and  temperature  distributions.  The  program  was 
entered  with  21  values  of  both  heat  flux  and  temperature  spaced  at 
equal  increments  of  r*  of  0.05.  Values  for  r*  increments  of  0.025 
were  obtained  by  linear  interpolation.  It  was  found  necessary  to 
limit  the  series  summation  to  the  first  ten  terms,  as  the  integration 
inaccuracy  became  unacceptable  for  the  higher  terms.  As  explained 
in  Section  V,  the  present  local  heat  transfer  results  are  inordinately 
sensitive  to  uncertainties  in  the  test  surface  temperature  distribution. 
For  this  reason,  the  program  was  not  modified  for  determination  of 
the  higher  order  terms. 

The  machine  program  has  two  procedures,  or  subsections.  The 
first  is  a  procedure  developed  for  determining  Jq(x)  ,  and  is  labeled 
JOX.  A  series  solution  is  used  for  x  less  than  10,  and  an  asympto¬ 
tic  expansion  for  x  greater  than  10.  The  accuracy  chosen  is  0.01 
per  cent  for  x  less  than  10  and  0.1  per  cent  for  x  greater  than  10. 
The  second  procedure  performs  the  integration  in  [  D  24],  and  is 
labeled  INT.  The  integration  scheme  is  a  simple  trapezoidal  inte¬ 
gration  . 

The  following  table  gives  the  translation  between  the  present 
nomenclature  and  that  used  in  the  computer  program. 
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A 

a  ,  the  zeros  of  J^(x) 

BRO 

b/r  ,  ratio  of  disk  thickness  to  radius 

COND 

disk  thermal  conductivity 

DR 

increments  of  r*  for  which  the  heat  flux  is 
to  be  determined 

DS 

increments  of  r*  for  performing  the  numeri¬ 
cal  integration 

ERR 

allowable  error  in  summing  the  series 

ING 

value  of  the  integral  in  [D  24  ] 

JO 

I0M 

JOA 

2 

Jo(am> 

M 

number  of  terms  summed 

N 

number  of  data  points  entered 

QB 

qb"(r*)(  heat  meter  section  heat  flux 

QO 

qo"(r*),  test  surface  heat  flux 

R 

r*  at  which  heat  flux  is  evaluated 

RO 

disk  radius 

RUN 

Run  number 

TEO 

t  (r*) ,  test  surface  temperature  distribution 

The  machine  program  follows. 
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COMMENT  EVACUATION  OF  QO  BY  SERIES* 

INTEGER  I  *J*K,M,N  .RUNS 

ARRAY  TIBI  50  >♦  TED  (50).  QBs'50)*  I  NG  (15)$ 

ARRAY  A  t 1 0  >  =  (0.0,  3.83171*  7.01559*  10.17347,  13.32369,  16.47063* 

19.61586.  22  *  76  008  *  25.90367.  29.04683)5 
ARRAY  JOA(IO)  »  (1.0*  0.162215*  0.0900695,  0.0623525*  0.0476808, 
.0378901,  0 © 0324228  »  0.0279507,  0.0245624*  0.0219081)* 

PROCEDURE  JOX(X$JO)$ 

BEGIN 

EITHER  IF  X  LSS  10. OS  GO  SML * 

OTHERWISE*  GO  LRGS 
SML..  F  »  X  «  X / 4$ 

SUMl  =*  TRM  a  R  =  1.0* 

UNTIL!  ABsrCRM /SUM  1  t  LSS  0.0001  \% 

BEGIN 

TRM  *  -TRM.E/ «R.R »* 

SUMl  *  SUMl  +  TRM* 

R  =  R  1.0$ 

END* 

JO  *  SUMl* 

RETURNS 

lRG.o  JO  *  0.79788H COSKX-O. 7854)  +  ( S I N ( X-0 . 7854 )  ) / ( 8X )  ) /SQR T < X ) * 
RETURN  End  JOXU’S 

PROCEDURE  INT!DS*N*Q( ),A**JOX<*» )* 

COMMENT  INTEGRATION  OF  (S  0  JO ( A  S))* 

BEGIN 

INTEGER  I*  N$ 

EITHER  IF  A  EOL  0.0* 

GO  ALPHO* 

OTHERWISE* 


GO  ALPHM* 


ALPHO..SUMO  a  0.5  (Q<N>>$ 

FOR  i  a*  C2»1»N“U$ 

sumo  *  sumo  +  { i-i  )DS<oi  i  ns 

I  NT  f  )  »  SUMO*  DSS 

RETURNS 

ALPHM..JOX( ASJO) $ 

SUMM  *  0.5 «Q( Nl >  JOS 
FOR  I  *  (2*1 »N-1 ) $ 

BEGIN 

ARG  -  (I-IHOS)AS 
JOX  (  A  RGS  JO  'i  $ 

SUMM  «  SUMM  -4  (  1-1  >DS(Q(  I  I)  JOS 
ENDS 

INT ( )  “  SUMM.DSS 

RETURN  END  INT( )$ 

COND  «  72. OS  BRO  *  0.2109$  R0  *  0.6649$ 

START ..READ( SSDATAjS 

INPUT  DATA(DS*N»DR.M*ERR*RUN»  FOR  I  «  ( 1 . 1 » N ) SQB l I) * 

FOR  I  *=  (  1  .1  *N)STEO(  I>  >$ 

TERM1  *  2.INT(DS»N  *QB ( ) *  0 • OSS JOX  (  $  )  )  $ 

WRITE  ($$ANSl»HDGl)$ 

OUTPUT  ANSKTERMDS 

FORMAT  HDG1 ( B5  »*J  =  i*.B24.*TERM  =  *  .  F 1 2 . 6 » W4 ) $ 

COMMENT  COMPUTATION  FOR  R  *  0$ 

SUMO  *  0.0$  TERMO  *  TERM1 $ 

I  =  1$ 

FOR  J  *  (2*1*10)$ 

BEGIN 

SUMO  *  SUMO  +  TERMOS 
FOR  K  *  (1*2*  2N  —  1 ) $ 

TIB  UO  a  QB (  ( K+l 5/25  -*■  A ( J ) . COND ( S I NH ( A ( J » . BRO ) I 


T£0  <  f  M  I  )  / 2  >  /ROS 


FOR  K  ®  <2*2*2N-2J$ 

TIB  IK)  *  <  T I  8  i  K"  I  t  4  TIBCK  +  I  >  ) /2.0$ 

INGOJ  «  INrfDS/2*0*2N-l»TIB(  )  *A{  J)SSJOX(S)  >* 

TERMO  *  2  f ING ( J  5 ) / (  fc1  JOA  (Ji  I  COSH  {  BRO«  A  (  J  )  J  )  $ 

WRITE  C  SSANS2  »HDG2  >  $ 

OUTPUT  ANS2  <  J  » I NG i J t »  T  ERMO ) $ 

FORMAT  HOG2 l 85  »* J  ♦ I  2  *B3  *  *  I NG  ***Fi3.6»B3. 

♦  TERM  a*  *  F13  »o  *  W4 ' S 

ENDS 

QO  *  SUMO  +  TERMO* 

WRITE  i SSAN53.HDG3 5 S 
OUTPUT  A  N  S  3  i QO  »  R  UN  I  $ 

FORMAT  HDG3<6lO.*R  *  0*,Bi6.*QO  »* , F 13 . 6 . B5 * *RUN* > I 4 » W4 i $ 
COMMENT  COMPUTATION  FOR  R  GREATER  THAN  OS 
FOR  I  »  ( 2 ♦ 1 »M ) $ 

BEGIN 

5UMM  «  0.0S  TERMM  *  TERMlS 
J  “2S 

ADD* • SUMM  =  SUMM  4  TERMMS 
IF  J  EQL  Xi$  GO  QUITS 
EITHER  IF  ABSf TERMM/SUMMi  GTR  ERRS 
BEGIN 

ARG  **  (Aij)l(W)  DRS 
JOX«  ARGSJ0  5S 

TERMM  «  2JO. ING< J) I JOAl J» )COSH(BRO.A( J) }) S 
WRITE  i SSANS4  *  HDG4 I  $ 

OUTPUT  ANS4I J*1NG( J» » TERMM )$ 

FORMAT  HDG4  t  B  5  »  * J  » * , 5 2 « B 3 ** I NG  **,F13.6*B3, 
*TERM  «*»F13.6*W4i S 


J  +  1$ 


60  ADOS 
END* 

OTHERWISE* 

QUIf. .00  fc  SUMM* 

WRITE  (  $*aNSS  »HDG5  If  * 

OUTPUT  ANS5I I I-t SDR»00*RUN »* 

FORMAT  HD65«B10»*R  • F8 ,1 » B 10 » *00  ®**E13«6* 

85  »*RUN*  » I  4  » W4 ) * 

END* 

60  STARTS 


FINISH* 


Appendix  Heat  Meter  Calibration 

The  heat  meters  used  on  the  rotating  disk  are  those  manufactured 
by  National  Instrument  Laboratories,  Inc.,  Washington,  D.C.  The 
heat  meter  consists  of  a  fine  mesh  copper  screen  bonded  to  either 
side  of  a  silver-tellerium  wafer.  Copper  leads  are  soldered  to  the 
screens .  The  temperature  difference  between  the  two  copper  and 
silver-tellerium  junctions  due  to  heat  flow  normal  to  the  wafer  pro¬ 
duces  a  net  EMF  output  proportional  to  the  heat  flux.  The  constant 
of  proportionality  must  be  determined  experimentally.  The  nominal 
heat  meter  specifications  are: 

Catalog  No.  HF  -  1  -  C 

Diameter  0.400" 


Thickness 

K  (Calibration  factor) 
Thermal  resistance 


The  heat  meters  were  calibrated  after  the  disk  was  bonded  to 
gether.  The  disk  was  inverted  and  placed  over  a  16"  diameter 
electrical  heater.  A  cooling  flow  of  water  was  directed  over  the 
test  surface  (the  upper  surface  during  calibration).  The  direction 
of  the  heat  flow  through  the  disk  was  then  the  same  for  both  cali¬ 
bration  and  testing.  The  disk  heat  flux  was  determined  from  the 
electrical  power  input  to  the  heater.  Energy  balance  tests 
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established  the  heat  loss  at  less  than  0.5  per  cent.  The  relatively 
thick  aluminum  disks  insured  a  uniform  heat  flux  through  the  heat 
meter  section  (the  protective  disk  was  0.375"  thick  during  calibra¬ 
tion  and  later  machined  to  0.065"). 

Model  tests  established  that  the  heat  meter  output  was  not 
dependent  on  the  temperature  level  of  both  junctions  over  the  cali¬ 
bration  range.  The  EMF  -  temperature  relation  is  only  slightly 
temperature  dependent,  and  a  single  junction  temperature  charac¬ 
terized  the  heat  meter  output.  The  maximum  operating  temperature, 
dictated  by  the  strength  of  the  bond,  is  170  F.  The  maximum  cali¬ 
bration  heat  flux  was  4000  BTU/hr  ft^  . 

The  heat  meter  calibration  factor  is  defined  as  the  heat  flux 
per  millivolt  of  electrical  output.  The  calibration  factors  for  the 
disk  heat  meters  are  presented  in  Fig.  E-l.  For  heat  meters  1,  3, 

4  and  5,  the  calibration  factor  varies  linearly  with  temperature. 

The  effect  of  temperature  is  slight;  a  100  F  rise  decreases  K  by 
about  8  per  cent. 

The  nonlinear  calibration  for  heat  meter  2  can  possibly  be 
explained  as  the  presence  of  a  bond  failure  between  the  heat  meter 
and  the  disk.  At  higher  temperatures  and/or  heat  fluxes,  a  gap 
between  the  heat  meter  and  the  disk  might  be  produced  by  differen¬ 
tial  thermal  expansion.  The  heat  flux  through  the  heat  meter,  and, 
consequently,  the  output  of  the  heat  meter,  would  be  reduced, 
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yielding  a  higher  value  of  K.  This  calibration  did  not  appear  to 
hold  during  testing,  as  the  calculated  heat  flux  was  about  twice 
that  measured  by  any  other  heat  meter.  An  adjusted  characteristic 
passing  through  the  two  lowest  temperature  points  was  used.  This 
gave  heat  fluxes  in  good  agreement  with  those  obtained  by  plotting 
the  other  heat  fluxes  as  a  function  of  radius  and  fairing  a  smooth 
curve  through  the  data . 


FIGURE  E  -  1 


HEAT  METER  CALIBRATION  FACTOR,  K,  AS  A 
FUNCTION  OF  HEAT  METER  TEMPERATURE 

(a)  Heat  Meters  No.  1,  3,  4,  and  5 


(b)  Heat  Meter  No.  2 


750 


(o) 


(b) 
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Appendix  F  Derivation  of  Integral  Equations 

The  integral  equations  for  the  fluid  dynamics  and  heat  transfer 

j 

from  a  rotating  disk  are  derived  from  the  governing  differential 
equations.  The  major  assumptions  are: 

1.  The  flow  and  heat  transfer  are  steady  and  axisymmetric . 

2.  The  fluid  pressure  is  a  function  of  radius  only. 

3;  The  fluid  behaves  as  an  uncompressed  gas. 

4.  The  flow  is  turbulent. 

5.  The  heat  transfer  is  in  the  axial  (z)  direction. 

The  governing  differential  equations  are: 

Continuity  equation: 

u  +  _u__  +  9  w  =  o  [FI] 

3  r  r  3  z 

Momentum  equation  in  the  radial  direction: 

„  2 

gc  3  fr  -  g^,  dp  =  u  3  u  +  w  3_u  v_  [F2] 

P  3  z  p  dr  3  r  3  z  r 


Energy  equation: 


-  1 

3  q"  = 

U  3  t 

+  w  3 

_t_ 

p  cp 

3  z 

3  r 

8 

z 

boundary  conditions 

are: 

At  z  =  0: 

u  =  0 , 

v  =r  w  , 

w  =  0 , 

Tr  =  Ttw 

t  =  *w 

f  Qn  =  Qw 

At  z  =  5: 

u  =  u6 

'  v  =  V6  ' 

w  =  w 

6  '  t  r  =  0 

At  z  =  d: 

t  =  t  A 

,  q”  =  0 
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[F3] 


l 


The  continuity  equation  [  F  1]  is  first  integrated  to  evaluate 
the  axial  velocity  component,  w,  at  the  edge  of  the  hydrodynamic 


boundary  layer,  <$ 


6 

^<5 

1 

[  f3  u 

+  u  ^ 

d  z 

+  i  3  w  dz 

— 

0 

J  h  r 

r  J 

J  3  z 

r\ 

wj  =  -  P_l_  5  ( r  u )  dz 

J  r  3  r 


Using  the  mathematical  relation  for  differentiation  under  the 


integral  sign,  [F  4]  can  be  written  as: 

r6 

W6  =  ~  J_  d_  r  u  dz 

r  dr 


u<s  d  6 
dr 


The  momentum  differential  equation  [  F  2  ]  is  now  integrated 
over  the  hydrodynamic  boundary  layer: 


gc  3  Tr  dz 

P  3  z 


gc  dz 

T-  dr 


-6  r~6 

u  5  u  dz  +  w  3  u  dz 
3  r  3  z 


v2  dz 


Considering  the  terms  separately  and  employing  the  boundary 
conditions: 


£c  3_Ic  dz 
o  P  3  z 

r  6 

-  gc  dp  dz 

JT  dr 


gc  T  r„ 


£c 

p  dr 


W  9u  dz 
32 


w  u).  -  w  u) 

O  O 


u  3w  dz 
3z 


Using  the  continuity  equation  [F  1]  and  f  5],  [  F  9  ]can  be 
written  as: 


w  3  u  dz 
a  z 


-  ufi  d  r  u  dz 
r  dr  ^ 


u  ^  d  6  + 


u  3  ( r  u  )  dz 
r  3  r 


t  F10  ] 


[F  6  ]  can  then  be  written  as: 


r6 

£ 

H  1 

+ 

6  dp 

= 

u  a  u 

+  u  a  (r  u) 

dr 

j 

L  ^7 

3  r 

u6  ru  dz  -  u6  d  6  -  v  dz 

r  dr  J  dr  J  "T 

°  o 


tFll] 


The  first  integral  on  the  right  hand  side  can  be  combined: 


u  3u  + 
"ar~ 


u  a  ( r  u  ) 

T  Fr  _ 


2 

dz  =  J.  3  (r  u  )  dz 
-  r  3  r 


[F12] 


Again  using  the  relation  for  differentiation  under  the  integral  sign: 


!  3 ( r  u  )  dz 
r  3  r 


[  2  2 
_1_  _d_  ru  dz  +  u6  d6 

r  dr]  dr 
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[F13] 


Employing  [F  13  ] ,  the  integral  momentum  equation  becomes: 


Normalizing  the  coordinates  and  <5  by  the  disk  radius  yields: 


The  energy  integral  equation  is  obtained  by  integrating  the 
energy  differential  equation  [  F  3  ]  over  the  thermal  boundary 
layer: 


Considering  the  terms  separately: 


[FIT] 


w  t)  A 


w  t)o 


dz 

[F18] 
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